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SUMMARY 

A  uniiicd  math emailcal  account  is  given  of  the  available  knowledge, 
together  with  additional  new  results  of  the  entropy  di st I'ibution  through 
a  normal  shock  wave  The  most  ’miable  leainre  ofth.-'s  d*strj button  is  tlie 
fact  that  as  long  as  heal  conductivity  is  present  the  entropy  will  first 
increase  within  the  shock  until  it  j'eaclies  a  maximuni  valiie  at  a  certain 
point  inside  ot  the  shock,  and  then  diminishes  to  its  Unal  v alu e  behind  the 
sliock.  A  physical  drscussaon  oi  the  results  is  given  in  add>tion  16  a  review 
of  the  phenomena  not  usually  uicliicU'd  in  the  analysid  oi  shviuk  wave  striiilure 
but  whieli  er-tcr  when  tiie  sti'ength  oi  the  shock  wave  :  s  suffici  eniJy  great  so 
that  tlie 'chemical  rear  Tons  lakf'  piacn 
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A  systematic  review  of  classical  shock  wave  structure  accoiding 
to  the  Navier-Stokes  equations  is  included  here  together  with  a  discussion 
of  the  physical  validity  of  these  equations.  The  structure  of.  and  the 
entropy  distribution  within  weak  shock  waves  in  general,  and  shock  waves 
of  arbitrary  strength  with  Prandtl  numbers  of  0,  3/4  and  co  are  analyzed 
in  detail,  together  with  qualitative  results  for  shock  waves  in  general, 

The  case  of  a  shock  wave  in  a  fluid  with  heat  conduction  but  without  viscosity, 
affords  an  example  of  a  system  within  which  a  discontinuous  change  of  slate 
to  a  lower  entropy  occurs. 
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I  -  INTRODUCTION 


For  roughly  tlie  past  twenty  years  tlie  distribution  of  the  entropy 
through  a  normal  shock  wave  has  been  of  interest  to  a  variety  of  investi¬ 
gators.  The  most  noteworthy  aspect  of  the  entropy  di.stribution  is  that 
although  the  velocity,  specific  volume,  temperature  and  pressure  of  a 
gas  all  vary  monotonically  ac ross' the  shock  structure,  tlie  entropy, 
instead  of  also  varying  (in  particular,  increasing)  monotonically  throughout 
the  shock,  increases  only  at  first,  passes  through  a  maximum  within  the 
shock,  and  then  actually  diminishes  to  its  final  value  behind  the  shock. 

This  phenomenon  has  proved  of  such  interest  that  a  perusal  of  the 
literature  indicates  that  this  .fact  appears  to  have  been  discovered 
independently  at  least  lour  dilferenl  times  in  the  period  from  1944  to 
196],  In  view'  of  this  situation  and  in  view'  of  the  inlrirusic  interest  in 
the  entropy  distribution,  the  purpo.se  of  this  paper  i.s  to  present  a 


systtnidlic  investigation  of  the  distribution  of  the  eiifrojiy  tliroughout  the 
the  structure  of  a  normal  shock  wave,  including  es.sentially  all  knowTedgc 
of  this  distribution  that  ha.s  been  heretofore  obtained,  in  addition  to  new' 
results  obtained  here  by  the  present  authors.  To  make  this  account  as 
.self-contained  as  possible,  a  unified  matliematicai  i.>rrs enlation  i.s  first 
given  of  shock  wave  structure.,  including  the  varioms  known  exact  clo.sccl- 


lonn  solutions  for  various  special  cases,  namely  ex- cemol  v  weak  shocks 
ill  gCncial  and  shock  wave.s  of  arbit  rary  strengtli  with  Prandti  mimhers 
of  0,  3/4,  and  co.  Because  of  the  important  role  played  by  tiie  heat 


C.onduc'iiviiy  in  the 
is  di .sc u ss I'd  here  i 


eiittopy  distribution,  Uic  ta.se  of  zero  Prandti  nutnbes’ 
a  .?omcw'1i.a!  greater  detail  Uiau  ordinarily  found  in  the 


literature.  The  entire  analysis  in  this  paper  is  based  on  the  Navier- 
Stokes' equations  of  flow;  accordingly,  a  review  is  given  here  of  the 
physical  validity  of  these,  equations  for  the  shock- structure  problem. 
Finally,  additional  physical  effects  such  as  dissociation  and  ionization, 
which  arise  in  strong  shock  waves,  are  briefly  discussed. 

Historically,  it  may  be  noted  that  Rayleigh^’  ^  in  1910  obtained 
exact  solutions  for  the  shock  structure  for  the  special  cases  of  a  gas 
without  heat  conductivity,  but  with  viscosity,  and  of  a  gas  without 

3,4 

viscosity,  but  with  heat  conductivity.  At  about  the  same  time,  Taylor  ’ 

obtained  a  solution  for  the  structure  of  weak  shock  waves,  including 

5 

both  viscosity  and  heat  conduction.  In  1922,  Becker  solved  the  shock- 

structure  problem  for  a  Prandtl  number  of  3/4;  this  is  still  the  only 

case  for  which  a  simple,  closed-form  exact  .solution  foi"  a  shock  wave 

of  arbitrary  Strength,  including  both  viscosity  and  heat  conduction, 

exists.  This  ca.se  has  also  been  considered  by  Roy^’  Thoma.s'‘, 

10,6 

Morduchow  and  Libby'',  and  Puckett  and  Stewart  ,  Hoy  and  Thomas 
considered  the  case  of  viscosity  cocfric.ient  U  proportional  to  the  square 
root  of  the  absolute  temperature  T,  while  Morduchow  and  Libby  con¬ 
sidered  the  more  general  ca.se  of  (a>0)  ReferenetJs  8  and  9,  in 

particular,  showed  the  important  influence  of  variable  viscosity  and 
heat  conductivity  in  thickening  the  shock  wave,.  A  ma thematicai  analysis 
of  shock-wave  structure  under  more  general  conditions  (e.  g.  ,  general 
Prandt!  number)  was  subscciuontly  made  by  von  Miscs  (reierence  i!)  and, 
i ndepicndcntl y ,  by  Gilbarg^^, 

Concerning  t!ir,  entropy  distribution  v/ilhin  t.hc  shock-wave 


structure,  it  appears  that  it  was  first  found  in  1944  by  Roy^,  considering 
the  case  of  a  Prandtl  number  of  3/4,  that  the  entropy  will  not  increase 
monotonically  across  a  shock,  but  will  pass  tli rough  a  maximum  at  a 

Q 

certain  point  within  the  shock  wave,  Morduchow  and  Libby'’  subsequently 
also  discovered  (independently  of  reference  6)  the  non-monotonic  behavior 
of  the  entropy  for  this  case,  and  showed  specific  curves  for  the  entropy 
distribution.  Zeldovich  (quoted  in  reference  13)  in  1946  found  a  similar 
behavior  of  the  entropy  for  a  sufficiently  weak  shock  wave  in  a  gas  wntli 
thermal  conductivity  but  without  viscosity,  and  Landau  and  Lifshitz 
(quoted  in  reference  13;  see  also  reference  14)  in  1954  noticed  a  similar 
behavior  for  extremely  weak  shock  waves  in  general. 

The  quite  recent  interest  in  the  entropy  distribution  through  a  shock 
is  evidenced  in  papers  by  Roy”^,  Ackeret'^,  Golitsyn  and  Staniukovich  ^ 
and  Serrin  and  Whang^^.  Serrin  and  Whang  (who  appear  to  have  been 


unaware  of  any  previous  finding.?  on  tlic  entropy)  prove  that  the  entropy 
will  vary  non-monotonically  across  a  shock  wave  under  quite  general 
conditions.  It  will  be  seen  here  that  the  most  important  condition  is  that 
the  gas  have  a  positive  thermal  conductivity  coefficient. 

The:  analysis  pertains  explicitly  to  normal'  shock  waves,  Howevor, 
it.  is  noted  that  by  adding  a  constant  velocity  normal  to  the  x  axi.s,  i.c,  ,  a 
V  vcdocily  sucli  shat  v=  constant,  the  structure  of  plane  oblique  shock  wavc.s 
is  implicitly  being  considered. 

9  lie  author.?  hereby  e.xpress  their  thanks  to  Herbert  Fox  for 

his  aid  in  the  calculations  and  to  Professor  M.H.  Bloom  for  hi,s  u.scfu] 
discussion.?. 


II  -  basic  equations  of  flow 


.1  ne  cttissrcai  equations  of  continuity,  momentum  and  energy  fo: 
the  onc-dimonsionai  .steady  H o>.v.  of  a  cotnpresaibl  c.,  viscous,  h.cat- 


conducting  fluid  (references  Zand  17)  are  respectively; 


(pu)'  =  0 

(1) 

puu  '  +  p  '  -  {4/3)(mu')  '  =  0 

(2) 

puE'  +  pu'  -  (4/3)].; u  -  (kT  ')  '  =  0 

(3) 

It  will  be  assumed  throughout  the  analysis  that  the  fluid  is  homogeneous ,  and 
a  peifect  gas,  i.  e.  ,  that  it  obeys  the  ideal  gas  law  and  has  constant 
specific  heats.  Later  in  connection  with  the  analysis  of  strong  shock 
waves  a  more  accurate  way  of  describing  equations  will  be  presented; 
several  of  these  assumptions  will  be  removed.  Thins,  here  llie  internal 
energy  E  is  assumed  to  be  given  by 

^  (4) 

while  the  equation  of  state  is  assumed  as 

P  =  pHT  (5) 

From  Eq,  (])  it  follows  that 

pu  =  m  ((,) 

wlicte  ni  is  a  constant.  Eq.  (o)  sliows  that  the  specific  volume,  v  {  =  l/p) 
wall  vary  exactly  like  the  flow  velocity  u,  .since  v  =  (l/?n)u.  The  entropy 
S  pe !  unit  mas.Sj  oniilting  an  additive  constant,  may  be,  WM'iUen  a.s 

5  =  log  T  -  R  log  0  (7) 


In  view  o£  Eq.  (6),  S  can  also  be  expressed  here  as 

S  -  log  T  +  R  log  u  (8) 

From  Eq.  (8),  in  conjunction  with  Eqs.  (4)  -  (6),  it  follow.?  readily  that 
puTS  ^  =  puE  ^  +  pu^.  Hence  the  energy  equation  (3)  can  also  be  written 
in  the  form 

OUTS'  =  {kT') '  +  (4/3)mu'^  (9) 

Using  (6),  Eq.  (2)  can  be  integrated  once  to  yield: 

p+  mu  -  (4/3 )mu'  =  (10) 

where  Is  a  cornstanl.  Substituting  for  p  according  to  Eq.  (lO)  into 
Eq.  (3).  it  is  found  that  the  latter  can  be  integrated  once  to  yield: 

mE  +  C^u  -  (m/2)u^  -  kT'  =  (11) 

where  C  is  a  constant.  Since 

3 

p  =  mRT/u  ( 1  2) 

Eqs.  ( 1  0)  and  (11)  can  be  written  in  ilie  .fo rin 

pu  '  =  lv!(n.  T)  (13) 

kT'  =  L(u,T)  (14) 

where  M  and  h  are  given  functions  of  n  and  T,  and  do  not  coniain  p  or 
k .  First  order  equations  ol  the  form  (13)  and  (14)  which  determine 
directly  T  vs.n.  have  been  u.sed  io  analy ve  quail iativel y  the  nature  of 


the  general  shock  wave  solutions  (references  11  and  1  Z)  , 


As  an  alternative  to  Eq?.  ( 1  3)  and  (l4)  another  useful  equation 
can  be  obtained  by  writing  T  =  {pu)/(mR)  and  using  the  expression  for  p 
in  terms  of  u  and  according  to  Eq,  (10),  Substituting  then  for  T  into 
Eq.  (11)  ,  the  following  differential  equation  for  the  velocity  u  vs.  x  is 
obtained 


3(Y-1) 


L.  K  ...  . 

_i_,  ]  -  — ^1_[u(hu  ')  '  *•  —  '  1 

Rm  3Rm  3 


(15) 


Y 

Y-l 


C  u  + 


m(Y+  1) 
Z(Y-l) 


Here  the  gas  relations  R  =  c  -c  and  y  =  c  /c  iuivc  been  used, 

py  [D '  V 

Finally,  another  fruitful  equation  can  be  obtained  by  noting  that 
according  to  Fq,  ( I  Z)  p^=  (mRT^/u)  -  (pu^/u),  s\ibsiiluling  for  p'  into 
Eq,  (Z)  and  solving  for  pu  \  Then  pvitting  this  expression  for  pu  ^  into 
Eq.  (3),  an  equation  is  obtained  which  can  be  integrated  once  (cf,  refer¬ 
ence  9),  to  yield 

(u  /Z)  I-  -  [k/mc^^][(4/3)Pr  ui/  +  ']  =  C  (lb) 

w'i'icre  C  is  a  constant  and  pr  is  the  Prandtl  nuinbcr  defined  Ity 

Pr-uc,,^/k  (17) 


Eqs,  (.3),  (9),  and  (16)  constitute 
equaiioius  Eqs,  (  3)  aiid  (9)  shv’W 


three  allcrnalive  forms  of  I'ne  energy 
th'.‘  role  of  heat  conductivity  ,md  viscous 


dissipation  in  the  determination  of  an  energy  balance.  The  importance 
of  Ea.  (16)  is  due  to  the  appearance  there  of  the  "total"  ( or"stagnation") 
enthalpy[(u^/2)  +  c^T]. 

All  of  the  preceding  equations  are  valid  for  arbitrary  and  variable 
viscosity  and  heat  conductivity  coefficients  IJ  and  k. 

Ill  -  SHOCK- WAVE  STRUCTURE 

A  unified  mathematical  account  of  the  solutions  lor  shock-wave 
structure,  based  on  the  set  of  equations  in  II,  will  be  given  here. 
Discontinuous  Shock  Wave:  1J  =  0,  k-Q 

If  viscosity  and  heal  c onduction  a  re  neglected,  i.  e,,  u  =  0  and 
k  =  0,  then  Eqs,  (6),  (10),  (16)  and  (5)  become  the  following  a1  gebraic 
equations 


pu  =  m 


p  ■)-  piE  -  C 


(18) 


p  =  CRT 


The.sc  equations  indicate  a  solution  of  the  form  p  =  con.sl,  ,  p?  const.  , 

I  -const.,  u=const  However,  if  one  solve.s  for  these  quantiiie.s  in 
term.s  of  tlie  consiants  m,  C_  and  it  is  tound  that  two  posKible  sets 
01  values  for  tlsesc  quantities  exist  The.  physical  interpretation  i  .s  lliai 

at  some  (here,  arbitrarv)  point  x  Z  there  can  be  a  djscontinuitv  surface 
across  wliirh  each  quanUtv  nimps  fro.m  one  constant  vaiue  to  the  other 


constant  value.  This  discontinuity  is  a  "normal  shock  wave/'  and  the 
solution  of  Eqs.  (18)  maybe  considered  to  be 


u  =  u  , 
1 

p  =  p, . 

P  =  P,  , 

T  ^  T 

\ 

(x  <  e) 

U  r  u  , 
2 

P=P3, 

P  =  P  , 
2 

T--T3 

(x  >  ?) 

(19) 


where  u  ,  u  ,  p  ,  d.,,  etc,  are  constants  (Fig.  1),  The  "strength"  of  the 

shock  wave  is  measured  by  the  magnitude  of  the  jump,  such  as  p  /p  , 

2  1 

The  relationship  between  conditions  ("l")  ahead  of  the  shock  and  con¬ 
ditions  ("2")  behind  it  can  be  found  by  noting  that  Eqs.  (18)  represent 
conservation  conditions  across  the  shock,  namely,  conservation  of 
flux  of  mass,  momentum,  and  energy,  respectively,  and  hence 


P  ri 


P 


p  H-  p  U  ^  p 

1  X  1  a 


+ 


3 

3 


c  T 
P 


c  •]• 
P 


(20) 


Eqs,  (20),  in  conjunction  with  ilie  gas  law,  can  be  used  to  solve  for  tin' 

ratios  u  /u_  ,  o_/p  ,  etc.  in  terms  of  a  sitigle  parameii-r  denoting  the 

sliock  strengtli.  One  convenient  parameter  is  the  Mach  number  M  su  /a 

1  Ml 

in  front  of  the  shock.  The  solution  of  Eqs,  (20)  (oihcr  tlian  the  "trivial" 
.solution  u  /u  -  1,  c-  /c  =li  etc  )  in  terms  of  M  (c.  g.  ,  reference  M)  is 

3  13  1  1  •  ^  ' 

//Cg  "  =  [(v-  DM  "  +  2]/[(y  -t  1)M^'] 


Pg/P.  =  f2vM,7(^  n]  -  [(v  -  })/(v  +  D]  (21) 

T^/T^  =  [2vM,;  -  (Y  -  l)K(v  -  DM./  +,  2i/f{v  t  1)"M.'M 


Eqs.  ( 2 !)•  c ons titute  one  form  of  the  well-known  "Rankine-Hugoniot 
relations”  between  quantities  in  front  of,  and  behind,  a  normal,  shock  wave. 

If  M  >  1  .  it  is  found  from  Eqs.  (21)  that  p_/p  ,  P„/P,  >  T  /T  >  1, 
while  1;  i-e.  ,  there  arises  a  compression  wave.  If  <  1,  the 

inequalities  are  reversed  (rarefaction  wave).  As  far  as  Eqs.  (20)  or 
(21)  are  concerned,  both  compression  and  rarefaction  waves  would  appear 
equally  possible.  The  entropy  S,  however,  also  jumps  across  the,  shock 
and  it  is  found  (e.  g.  ,  reference  14)  that  for  compression  waves, 

but  <  S  for  rarefaction  waves.  Consequently,  the  second  law  of 
thermodynamics  requiring  ,  may  be  invoked  here  to  show  that 

only  compression  waves  may  exist  physically,  and  that  a  normal  .shock 
can  therefore  occur  only  when  the  flow  Is  supersonic  (M^  >  1).  The  use 
of  the  second  law  of  tliermodynamic s  here  is  in  contra.st  with  the  results 
when  viscosity  and  heat  conductivity  are  taken  into  account.  In  that 
case,  as  will  be  subsequently  seen,  a  shock  -  st  ruclu re  solution  is 
mathematically  possible  only  for  >  1,  and  thus  the  second  law 
becomes  automatically  satisfied. 

Solutions  with  Viscosity  and  Heal  Conduction  in  General.  Boundary 
C  onditions  ■ 

When  ij;=k  =  0,  it  has  been  .seen  that  the  shock  "structure”  con- 
.sists  mo.reiy  of  a  discontinuity  surface  of  zero  thickness,.  In  the 
pre.scnce  of  visco-sity  and  heat  conduction,  however,  the  structure  may 
now  be  c.xpcct’ed  to  consi.stof  a  continuous  transition  from  state  1  to 
slate  2,  with  mo.st  of  the  change  occurring  in  a  small,  but  non-zero 
width  about  x  =  In  particular,  referring  to  Fig.  1,  it  may  he  expected 


0 


that  no»  .ha  actual  dis.ribu.iou  of  the  velocity  (say)  versus  the  d.stauce 
X  should  be  a  con.iuuous  curve  which  is  asymptotic  to  the  s.ra.ght  li„e 
U=u^  as  X  -  -CD  (corresponding  to  state  1)  and  asvmnto^ic  -  .  . 

liue  ueu^  as  x  -too  (correspondiug  to  slate  Z).  Cousequently ,  the 
appropriate  boundary  coudilious  to  be  .mposed  in  solving  the  differential 
equations  for  the  shock  structure  are  lhal  the  flows  be  uniform  a,  x  a  ±  co, 

..  e.  .  me  x-der.„ative  of  each  of  the  flow  variables  must  approach  aero  as 
X  V  ±®-  Moreover,  the  Ranki „e-Hugoniol  relations  will  hold  belween 
states  at  x  ^  _  oj.  These  relations  are  independent  of  the  viscosity 
and  heat  conducdv.ly  coefficients.*  Mises  and  (mdependentiy)  GilbarglM 
have  shown  mathematically  that  there  exists  one  and  only  „„e  sohm,,,,” 
namely  the  shock -structure  solution,  of  the  differenlial  equations  in  11 
satisfying  these  nnitorm  cond.lton,  at  x-h  on,  for  real  gases  ,n  genor.nl, 
"■i>h  posittvc  Vtscosity  and  hea,  conducvity  coeff.c cn.s,»  ,Mso,  cibarg'^ 

'  . .  **  R  h  both  independently  .approach 

zero  for  o  real  fluid,  the  KheSe’--  Ilvcm-  aonrnt,-ii,  i  r  ~~ 

^  .  1  e.  \  c.  1  a  pp  1  od h  e  s  i  h c  d i  ,s-  c  o  n  t i  r.u  ou  ,s  ( y  p e 

01  s-olution  (Fig.  !)  In,-  .k  =0.  AUhoiigh  the  equation.s  for  the  general 

case  must  be  solved  nnnterically  or  approxi rnately ,  .n  nombe,- 

mauVb^ltaKnmlnrm'o'mld'c^  ■ 

ditions  L  xs-o,, "^hcA  at;;  p’.;;'  ■  "’>■  '“i  "m  from 

thoHe  equ,aSions  .  '  Thi^  ^  eads  lo^Eor  Uor' W  i ’ "  ^  i 

that  at  points  wliere  the  flow  ‘mrif  m  V .  if  .should  be  nolecj 

conduction  t<nuns  tend  viscosity  and  heat- 

h?/d^;t“.AtLtu’qVM;!,°tn‘‘''  yf”"*  ">e  type  ( ,  .-,„d  |  ,,i,  whence 

a.  10.  .  in  the  iu,’h  ‘I'anV,'  The mtf 

and  satisfying  ihe  Ranid  ne  -  Hutson' o'  ’tH-  to.Krio^, 

this  difforeniial  eauation  '  j  e  ""  th  ^ -stncular  points  of 

Uu,T}  =  0  and  Mfu/l'lt  Thi  ;«,T)plane  where  both 

KOiutior,  .icnmng  ;lu-.se  (wo  pointpi  in'iHe'f!!'’ ‘rV  ^  unique 
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instructive  exact  solutions  have  been  obtained  for  certain  special  cases, 
and  these  will  now  be  given, 


Prandtl  Number  of  3/4 


If  Pi  has  a  constant  value  of  3/4,  then  Eq.  (16)  becomes  a  first- 


oider  linear  differential  equation  in  the  stagnation  enthalpy  [(u^/2}  +  c  TJ. 

The  solution  of  this  equation  for  uniform  conditions  at  x=+cd  (cf.  reference 
9)  is 


(uV2)h-  CpT  (22) 

w'liere  is  a  constant.  Thus,  (he  stagnation  enthalpy,  which  must  in  all 
cases  [see  Eqs,  (20)]  remain  the  same  behind  the  shock  (x  =  4-  cti)  ns  in 
front  of  the  shock  (x  =  mo),  will  in  this  case,  and  only  in  this  case,  also 
1  emain  constant  Ih  roughout  the  shock  structure.  Substi  Luting  i  nto 
Eq,  (10)  for  p  in  terms  of  u  according  to  Eqs.  (12)  and  (22).  and  solving 
for  pu  ,  one  finds  [ju'=F(u)/u,  where  F(u)  is  a  quadratic  in  u.  Since 
u  '=0  at  u  =u^  and  at  u  it  follows  that  F(u)  ~(u-u  )[u-u  ).  Thus 

it  is  found  that 


du  _  3(  y-l-  1 )  (u-ui  )(u-U3) 


it  o.’.y  be  remn  rked  Imre  that  if  the  conditions  of  uniformity  at  x  =  -m 
v-imuv'^''/f?  mathematically  a  consuiernble 

m  over  the  fuli  domain 

theT^r;  onm  b  ^-ererences  9  and  18),  The  present  authors  si, owed 

n'm.  !  ,  ■'“’■'-''’Ons  ,n  reference  9  but  did  not  discus.s  their 

Dobued  outTharfr"'’'l  Pi-ivate  communication 

at  1  n  im  V  T,  correspond  to  a  .source  or  sink  of  energy 

Covild  be  exliibiting  this  behavior 

<=0  ns  m  r  ^  restrict  excessively  the  gas  How 

of  eu..v,v  , cm addition 
.V  ...c  ^a-T,  u\  upstream  conduction. 


Eq.  (23)  holds  for  variable,  as  well  as  constant,  M  By  differentiali iig 
(iroplicitly)  the  right  side  of  Eq.  (23)  with  respect  to  a,  and  noting  that 

d/clx  -  [o/clu) ,  (.1  j  it  is  reacdily  found  that 


(tm  '■)  ' 


u  =  V  U ,  u 
1  2 


=  0 


(Z4) 


In  the  solutions  to  be  obtained  for  Pr.=  3/4,  the  origin  for  x  will  be  chosen 
at  the  point  of  inflection  of  the  vel  oc i ty- di  stanc  c  curve  For  constant  U. 

For  constant  p,  tlie  solution  of  Eq  {Z}) 


this  means  that  u  =  ViTV  at  x  =  0 

12 

is  then 


X  = 


8y  m  ry  '•‘l  .  .  -U 


3(y+  1)  m  Uj  -u 


[( - log  (- 


u,  -V  u,  U 
1  1  3 


,  r  *-‘3  .  .  .  tl-u„ 

-)  -  ( - )  log  ( - ^ - )] 


u  -u, 

I 


v^l  u  -u 

12  3 


(25) 


For  M  varying  in  a  presc  ribed  fashion  (e,  g.  ,  with  tempo  ralnre  T,  and 
hence  will,  u),  Eq.  (23)  can  be  solved  by  numerical  integration. 

From  £(}.  (23)  it  follows  that  a  continuous  solution  for  n(x) 
saiisCying  the  conditions  u  a t  x  =  -  m  and  u  =  at  x  .  ™  ,s  possible 

only  If  the  flow  .s  supersonic  In  the  region  ("1")  ahead  of  the  shock  One 
way  of  seeing  this  is  to  note  that  if  u  were  to  vary  continuously  from  a 
value  lo  a-hi^  value  there  must  be  some  region  where  n  is 
increasing  while  taking  on  some  values  between  u  and  u  .  ilowcwer 

1  9.  ■  ' 

tins  U  impccc.blc,  since  (.w  <  0  whenever  e  b„,wcen 

u..  and  hence  ,iceord,„g  lo  Eq  (2i|  „  0  whenever  „  us  belween  n  and 

\ 

Conricquenily,  musi  be  U.^s  than  u^  ,  and  this,  .nccording  to  the 
Rankine-Hugomot  relations  (21),  .mplms  M  >]  (cf  also  reference  O), 
Inns,  a  continuous  shock- structure  solution  is  maihcmaneniiy  possible 


only  for  a  snpersonrc  flow  in  front  of  the  ehock.  This  oonclnsion,  whioh 

has  been  establishee  here  (or  a  Prandtl  nnntber  of  3/4,  has  been  shown 

mathematically  to  hold  lor  all  positive  PrandH 

io-.-Uo..  .rUnrroeaa  (lelerence  il). 

Typical  curves  showing  the  distribution  of  vel  ocity  n /„,  (or 

spccific  volumel  and  of  temperature  T/T^  are  shown  in  Fig.  2, 

General  Prandtl  Number.  Pr-po. 

From  Fig.  2  il  is  seen  that  although  the  entire  transition  from 
state  1  to  state  2  occurs  over  an  infinite  domain  -co<it<cD,  by  far  the 
'nngest  portion  of  this  transition  occnrs  over  a  relatrvely  .small  distance 
-lied  the  shocit-wave  IhichnessS.  around  the  neen.or"  ,.vsb,  of  the  sboeh. 
Also,  the  velocity,  temperature.  Pressure  and  density  all  va ry  mono- 

tomcaUy  throughout  the  shoch,  the  velocity  decreasing  and  the  others 
■ncreasmg.  These  results  are  typical  of  the  shock- structure  .solutions 
general,  e.  g,  .  tor  general  Prandtl  number."  For  weak  shock  waves 
^»r  ouample,  it  will  bo  seen  below  tha,  the  form  of  the  shock- struct,, re 
•solution  remains  the  sanie  for  all  Pr,  ,he  latter,  tor  ,s  given  u, 
affect, „g  on,,  thickness  of  the  shock  wave,  M.ses'*  a„d  cnbarg'^ 
dare  shown  that  for  tl.e  shook- si  rue, „  re  soiuiion  in  general,  M(u,T)<0 
rind  L(u/r)>0  Ihrouohout  the  shorl-  Mr 

ui  UK.  shock,.  Hence,  irom  Eq.s.  [13}  and  (PI) 

itlollowfs  that  u  will  in  ycncral  decrease,  and  T  vil!  •  i  ■ 

C,  ana  t  will  m  general  increase 

monolonically  across  the  shock.  FromPo  c  ,, 

t-'q.  (0)  It,  then  lollov.'.-,  that  p  will, 

Vanous  speci.fic.  ciofinition.s  (cf  references  4  o 
ness  t  have  been  pro;ioscci  and  reed  t)  '  ’  ^  ^  iWck- 

<Mul-ua)/ldu^x(  most  common  is 

'  max 

ihc  Only  exception  is  the  case  of  i-tQ  k  y'O  fP,-  Mr, 

-■Giose  strength  exceeds  a  certain  anuw nt  n  ‘  ^  ^ 

subsequently,  .'’''ill  be  di  senssed 


inca-easemonotonically,  and  hence,  according  to  Eq,  (5),  the  pressure 
will  also  increase  monotonically  across  the  shock.  Thus  (exceptTor 
)  the  solution  foi  Pi  =3/4  may  be  regarded  as  the  prototype  of  the 
general  shock  solutions.’  Intact,  for  Prandtl  numbers  above  3/4,  the 
shock-wave  structure  changes  relatively  little.  This  can  be  seen 
explicitly  by  obtaining  the  shock- structu re  solution  for  the  extreme  case 

of  a  viscous  fluid  without  heat  conduction,  i  e.  ,  k  =  0,  p/O,  Pr  =  co. 

First,  It  should  be  noted  that  since  u'’=u"=:0  where  u  =  u  and 
the  right  side  of  Eq.  (15)  can  be  written  as 


C 

3 


C  u  + 

Y-J  '  2(Y-1) 


-  1) 
2(y-I) 


(u-U  ){u-U  } 
1  2 


(26) 


Consequently  Cor  k  =  0  Eq.  (15)  reduces  to 


dx  8  11 


Eq.  (Z7)  IS  the  same  as  the  corresponding  equation  [Eq.  (23)J  for 
P-3/d,  cseep.  that  the  right  side  of  (27)  no  longer  contains  a  f.ncor  of 
y  III  the  denoniinator.  Th,,s  implies  ,h.nt  ,f  t.  a  F(e)  i.s  ,he  solnlion  lor  the 
telocity  vs.  X  lor  Pr  =  3/4.  then  Ihe  solution  now  for  the  velocity  is 
d  =  F|yx).  Thus,  as  Prvarresfron,  3/4  loco,  the  »hock-wavo  thickness 
is  diminished  by  a  factor  of  1 /v  Shapiro  and  Kl.ne^^  h.ave  invesiigaied 
appro,x.ma,ely  ihe  effect  of  Prandt!  ntimber,  in  general, on  ,l,e  ihickneas 
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of  shock  waves  and  have  found  that  as  the  Prandtl  number  increases  the 
thickness  decreases,  this  effect  being  greatest  at  low  Prandtl  numbers. 
Weak  Shock  Waves. 

A  weak  shock  wave  means  here  one  in  which  [\i  -u^)  is  sufficiently 

small  so  that  only  the  lowest  necessary  powers  (namely,  second)  in 

(u^-u^)  need  be  retained  in  solving  the  equations.  Thus,  the  solution 

obtained  heie  will  be  valid  for  Mach  numbers  M  only  slightly  greater 

than  one.  This  solution  is  the  only  one  among  the  known  exact*^  closed - 

form  solutions  in  which  the  Prandtl  number  is  arbitrai-y.  Moreover,  as 

will  be  indicated  more  fully  in  the  next  section,  there  appears  little  doubt 

about  its  actual  physical  validity. 

To  the  degree  of  approximation  for  the  weak  shock  waves,  the 

va  nation  of  U  and  k  (c.  g.  with  T)  may  be  neglected.  Since  u  varies 

between  u^  and  ,  it  is  seen  from  Eq.  (26)  that  the  right  side  of  Eq.  (15) 

will  be  of  order  (u^-u^)".  Hence,  in  Eq.  (15)  u'  will  be  of  order 

'  rhciefore,  tlic  u  "  term  may  be  neglected.  Aforeov'cr, 

u  =(diP/du)u',  .  .  n"  is  of  order  u^^/{u  -u  )  or  of  order  (u  -u  )^  and 

1  3  13 

hence  the  u  term  may  also  be  neglected  in  Eci.  (15).  Finally,  to  second 
powers  in  [u  -u  ),  u  rnav  be  replaced  by  u  in  the  first  factor  of  u ''  in 
Eq,  (!5).  Thins,  Eq.  (15)  red\  ices  to  the  form 

U-'  =  (l/D)(u-u  )(u-uj  (28) 

°Of  course,  strictly  speaking,  the  solution  to  be  given  here  for  weak 
■shock.s  is  not  "exact,"  except  for  (u^ -u^)  -  0.  However,  for  (>,y  -u  ) 
sufficiently  small  it  will  be  a  rbit  ra  r\ly  c!  o.s  e  to  the  exact  solution.''^ 


where  D  is  a  positive  constant,  namely^ 


.1)] 

0  f v  +  n  t f  V- 1 )  p  ,, 

Fiom  Eq,  (Z8)  It  follows  thatu"=0  where  u  ^  (u^ +u^)/2.  Choosing  the 
origin  for  x  where  u"=0,  the  solution  of  Eq.  (28)  can  be  written  as 

o  -u 

^'-%-(~-Mtanh[(-^)x]  (30) 

2  2D 

whole  u^-(u)^^^q=(u^  +u^)/2,  The  temperature,  pre.ssuro,  and  density 
behave  quite  similarly  to  the  velocity,  since  by  Taylor's  scncs  the 
lemperaturc  (for  example)  will  be  T-T^^=  K(u-u^)  !  ,  who  re  U.(dTAiu)^^, 

andjt/0.  1'"  Thus,  to  the  degree  of  approximation  here.  T  vanes  linearly 
with  u.  Substituting  for  (u-u  J  and  evaluating  and  ;  so  (I, at  T  T  at 
^  -  m  and  T  =  at  x  =  +  ay  one  finds 

T^+T^  T^-T  u  -u 

-  =  ( - )  tanh  [(-i— i)  X  1  /an 


with  exactly  similar  results  for  the  density  and  the  pressure.  For  weak 
shock  waves  U  is  seen  (hat  the  form  of  (ho  shock  structure  renwunH  the 
same  for  all  u  and  k,  The  latter  affeei  only  the  constant  D  and  hmec, 

/.I  _  '  i- •  I  u\  Since  nk/a-.  =a,  .  and  hi  mu sf  he  r > ,v  t  i- 

p,  V,„,'l  =  I/lvM/)  ,nay  bf  D  bf'f/; 

Ko.'  ibP  oPivopy,  h&,. ,;=(dS/b„)_^=0  (Pf  Souio.,  IV  l,a„w) 


only  the  thickness  t  (see  footnote  5)  which  for  the 


weak  shock  waves  will 


t  =  4D/(u^-uJ  (3  2) 

Eqs.  (12)  and  (29)  show  that  an  increase  in  u  or  k  (or  both)  will  increase 
the  shock  thickness.  Eq.  (32)  shows  that  for  shock  waves  the 

thickness  decreases  rapidly  as  the  shock  strength  increases. 
of  the  Shock-Structure  .Sol, i, inn. 

According  to  the  shock- structu re  solutions  obtained  here,  the 
shock-wave  thickness,  except  tor  very  weak  shock  strengths,  will  bo 
on  the  order  of  themetn  free  path  of  the  gas  molecules  upstream  of  the 
This  can  be  seen  from  the  fact  that  the  natural  unit  6  (say)  for  x 
I's'ofcf.  e.g.,  Eq,  (25)]  i  s  6  =  M, /( n^ )  =  (tt^p^  a  ,  with  >  1 

For  a  standard  aeronauttcal  atmosph erc.i  at  om  atmosphere  prcsltre, 
6-(l.  69  X  10  /UJ  inches.  This  distance  is  nl  the  order  of  inagniliide 
of  the  mean  free  path  I  of  the  gas  molecules  the  kinetic  theory  of 
gases.  This  can  be  seen  not  only  numerically,  but  also  from  the  relation 
(e.g.,  reference  23)  .1;  =  ape  t,  where  J  =  average  random  (ihcrmal)  speed 
of  the  molecules  =  a,  and  a  is  a  constant  with  a  value  of  roughly  1/2, 

Thus  6-  f/|2MJ,  Specific  Ihcorclical  values  of  the  shock-wave 
thickness  can  be  found,  for  example,  in  references  9  and  22.  Becaua,, 
of  lhe.se  very  low  thicknesses,  Becker^,  who  assumed  conslatu  vl,,c„,,Uy 

_  .Vi9.7  1b/i„=,  ^ 

“V 1  1  4( !  TTsKTiD  I  Ji8  ft/sec,  '  ’ 


and  heaC;-c onducting  coefficients  and  found  thicknesses  even  lower  than 

a  mean  free  path  ahead  of  the  shock  for  moderate  Mach  numbers  M  , 

1 

originally  questioned  the  physical  validity  of  the  solutions  on  the  ground 
that  the  equations  used  here  are  valid  only  when  characteristic  distances 
along  the  flow  at  least  exceed  a  mean  free  path.  Subsequently,  Thomas® 
and  Morduchow  and  Libby^  showed  that  taking  into  account  the  actual 
increase  of  p  and  k  with  temperature  leads  to  significantly  higher 
thickne  s  se  s . 


The  question  of  the  physical  validity  of  the  basic  equations  used 
here  and  their  replacement  by  more  appropriate  equations  when  they  are 
not  valid,  is  doubtless  the  most  dit'ficult  in  the  problem  of  shock-wave 
stuictinc,  (Indeed,  if  the  basic  equations  were  known  to  be  physically 
completely  valid,  the  problem  of  the  classical  shock-wave  stnu  lure, 


i.e,,  without  effects  .such  as  di  s.sociatioii  or  iouizatioii,  could  be  con¬ 
sidered  as  essentially  solved,}  The  basic  equations,  namely,  the  Naviei- 
Stokes  equations,  used  here  are  the  continuum  equations  of  flow  of  a  gas, 
i.c.  ,  they  can  be  obtained  by  considering  the  gas  as  a  continuous  fluid, 
in  which  a  small  isolated  element  of  volume  i.s  to  be  in  equilibrium  under 
the  macroscopic  forces,  namely,  pressure,  shear  and  inertia  forces 
acting  on  this  element.  Moreover,  the  energy  equation  is  also  governed 
b\  mac  loscopic  quamit  i  c  s  ,ac  ti  on  a  volume  edemeut.  namelv,  the  inuu 
concluclecl  across  the  element  (bringing  in  the  temperature)  and  the  worK 
done  by  the  pressure  and  shear  forces.  The  Navier-Stokes  equations, 
in  particular,  also  involve  the  assumption  that  the  stress  components 
are  hnearly  related  to  the  rate-of-st  rain  compoueut  s ,  and  the  assumption 


that  the  hydrostatic  pressure  is  the  arithmetic  mean  of  the  three  normal 
sttesses,  makmg  the  bulk  (rate  of  dilation)  viscosity  coefficient  sero 
(  f.  reference  2).  The  continuum  equations  are  generally  considered  to 
lead  to  satisfactory  descriptions  of  the  flow  at  least  whenever  significant 
changes  of  state  occur  over  distances  L  considerably  greater  than  the 
mean  free  paths  of  the  gas  molecules.  Thus,  one  criterion  is  that  the 
Knudsen  number  (L/t)  must  be  at  least  of  a  certain  magnitude.  A 
c  atcd  cntoiion  for  the  validity  of  the  continuum  equations  has  been 
Ihe  number  of  collisions  experienced  by  a  molecule  during  its  transit 
across  the  shock;  Puckett  and  Stewart'“  have  concluded,  on  this  basis, 
lhat  in  the  shock  wave  there  are  sufficient  collisions  for  the  molecules 
10  .attain  equilibuum  for  their  translational  and  rotational  degrees  of 
feedom.  and  lhat  henoe  the  continuum  equations  should  be  approximately 
Such  t  .ileiia, which  are  relatively  simple  and  easy  to  apply,  may 
be  useful  in  the  absence  of  further  theoretical  and  experimental  infornio- 
lion  but  must  be  regarded  as  still  somewhat  arbtltary.  Also,  such 
crilona  usually  depend  strongly  on  ihe  poinl  in  ,|,e  flow  „ 
mean  fice  palh  is  calculated  |e.  g  ,  references  10  and  19).  The 
following  appear  at  present  to  be  the  oniy  true  criteria  for  delermiuing 
■ho  validity  of  the  continuum  equations  for  the  shock mvave  structure; 

(a)  How  do  Ihc-coutiiuiiim  sol„u„„a  compare  with  reliable  experi- 
mental  determinations  of  the  shock  .st  ructu  re  ? 

(b)  How  do  the  con.inuum  solutions  compare  with  siiffici  en.ly 

occura.c  solutions  of  the  Boitamann  cquahon  for  the  shock  -  st  ructurc 

problern? 


Both  criteria  (a)  and  (b)  are  ut  present  difficult  to  apply,  but  some 

progress  has  been  made  along  each  of!  these  directions.  Experimentally, 

limited  data  has  been  obtained  from  different  techniques  by  Hornig,  Greene 
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and  Cowan,  (references  20,  24-  26),  Sherman  and  Talbot  '  \  and  Ballard 
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and  Venable  .  In  discussing  the  experimental  results,  it  will  be  most 
suitable  to  review  first  the  theoretical  investigations  according  to  the 
kinetic  theory  of  gases  and  the  role  of  the  Navie r-Stokes  equations  in  this 
theory. 

Although  the  physical  validity  of  even  the  Boltzmann  equation  lias 
at  times  been  questioned  in  connection  with  shock  structure,  it  may  be 
assumed  at  present  that  this  equation  is  physically  valid  for  this  case, 
and  that  it  holds,  in  fact,  whenever  the  state  of  a  gas  does  not  change 
appreciably  over  a  distance  comparable  to  the  molecular  cl i a m e t c r 
(e,g.,  references  29"31)  .  In  the  Boltzmann  integro-'diffcrential  equation 
the  unknown  is  a  distribution  function  indicating  the  distribution  of 
velocities  among  the  gas  molecules  in  space  and  lime,  due  to  molecular 
collisions  and  into r- mol ecula r  forces,.  In  a  steady  equilibrium  state, 
i.  e.  ,  with  comparatively  small  or  negligible  gradients,  the  distribution 
function  is  of  the  -well-known  Maxwellian  type.  From  the  point  of  view 
of  kinetic  theory,  the  Navi  e  r-Stok  e  s  equations  can  be  regarded  as 
follows.  'riiC  Boltzmann  equation  can  be  multiplied  by  a  function  "(*?) 
of  the  velocity  vector  E  and  integrated  over  the  velocities.  Tlie.se  yield 
various  "moment"  equation^  In  particular,  if  -oiE)  i  .s  taken  as  1,  f, 
and  respectively,  and  average  (macroscopic)  quantities  .such  as  mean 
flow  velocity  u  and  density  p  0. re,  a ppropriat ely  defined  (e.  g,  rr-ference  30, 


31  and  19)  then  for  one-dimensional  steady  flow,  equations  of  the 
following  form  are  obtained 

(pu)'  =  0 

puu  '  +  (p+  t)  0  (33) 

pu(E  ''  +  UU  ')  +  [u(p+  r)] '  f  q  '=  0 

w'here  p  is  the  (hydrostatic)  pressure,  t  the  x-component  of  (shear)  stress, 
and  q  the  x-component  of  heat  flow.  It  is  easily  verified  that  the  set  of 
equations  (33)  is  exactly  equivalent  to  the  set  of  mass,  momentum,  and 
energy  equations  (1)  -  (3),  if  one  sets 

T  =  -(4/3)Mu'  ,  q  =  -kT'  (34) 

Thus,  the  Navier-Stokes  equations  can  be  regarded  as  arising  from 
Boltzmann's  equation,  provided  that  the  expressions  (34)  are  assumed 
for  the  stress  component  and  heat  flow,  i.e.  ,  t  (liiK-.arly),  proportional 
to  velocity  gradient  u'’,  and  q  proportional  to  temperature  gradient  T^. 
Indeed,  in  the  Enskog-Cha  pman  scheme  of  solving  th  c  B  o!  tzmann 
equation,  a  series  solution  in  powers  of  a  small,  pararncte.r  t  (the  m'Mn 
free  path.)  is  sought;  the  results  given  sequence  of  approximations  to  t 
and  q,and  the  first  approximation  for  t  and  q  is  tiiat  given  by  Eqs.  (34) 
(retercmcc  32).  Thus.,  the  Navi  er-Stoke.s  equations  arc  a  fir.'^l  approxi¬ 
mation  in  this  scheme.  It  should  be  noted  that  the  convergence  of  thi.s 
.scheme  in  general  In", .s  not  been  e.stabl i shed ,  and  in  fact  now  appoar.s 


r 
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even  doubtful,  Talbot  and  Sherman^  .  in  fact,  have  obtained  numerically 
a  solution  for  the  shock  structure  according  to  the  next  higher  approxima¬ 
tion  in  this  scheme,  namely,  the  Burnett  equations  and  have  found  (for 
Mach  numbers  up  to  2  roughly)  that  the  solution's  do  not  differ  signifi- 
cantly  from  the  Navier-Stokes  solutions  and  do  not  agree  any  better  with 
experimental  data.  Zoller^"^  had  previously  solved  equations  similar  to, 
but  not  identical  with,  the  Burnett  equations,  but  could  not  obtain  resvilts 
above  =  2,  36.  These  solutions  do  not  agree  as  well  with  experimental 
data  a.s  the  Navier-Stokes  or  Burnett  solutions  (reference  27),  An 

iterative  method  of  solution  of  the  Boltzmann  equation  has  been  con- 
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structed  by  Ikenberry  and  Truesdell  in  which,  again,  the  first  ite'i'ates 
yield  the  Navier-Stokes  equations.  Grad  (references  30  and  19)  has 
obtained  solutions  for  die  shock  structure  by  means  of  a  thirteen- 


moment  approximation  to  the  Boltzmann  equation.  The  solution  breaks 
down  for  a  Mach  number  of  1.  65  or  higher.  Wang  Chang^^  ha.s  obtained 
a  kinetic  theoiy  solution,  by  series,  lor  a  weak  .shock  wave.  The  con¬ 
vergence,  liowevcr,  is  slow  and  the  solution  breaks  down  for  shock 
uaves  ol  moderate  strength.  Mott -Sm.i lli ^ has  obtained  an  apj3roxi- 

male  solution  for  strong  shock  wave.s  by  assuming  the  distribution 
function  within  the  shock  to  be  a  (variable)  linear  combination  of  the 


Maxwellnan  ch  s!  ribuu  ons  in  Iron*,  ol  and  beliinc!  the  shock  This  solution 


'^Ikenberry  and  Tniesdeli^^  have  shown  that  in  tlic 


. n  Uie  special  case  of  a 

.simple  shearing  flow,  the  Navier-Stokes  solution  may  be  considered  a.« 
an  asymptotic  solution  and  cverv 
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higlicr  a pproxi ma ti oh ' j s  wor.sc  than  it 


Gi  ad  has  indit.ateci  that  the  convergence  of  this  solntion  may  be 
improved  by  expanding  in  a  different  parameter. 


T 


has  not  agreed  well  with  experiments  for  Mach  numbers  below  2^^’ 

however,  a  recent  modification  by  Muckenfuss^”^  to  take  into  account 

more  realistic  inte rmol ecula r  forces  has  led  to  improved  agreement 

with  experiment.  In  particular,  for  up  to  about  1.8  the  agreement 

with  experiment,  is  approximately  as  good  as  the  Navi er-Stokes  solution. 

but  at  =  2.  4  the  agreement  with  experiment  appears  appreciably  better 

for  the  kinetic-theory  solution.  A  further  modification  of  the  Mott-Smith 

type  of  solution  has  been  recently  made  by  Glansdorff^^  and  Ziei-ing, 

■tr*  -  39 

Ek  and  Koch  .  It  should  be  noted  that  for  very  weak  shocks,  the 
Navier-Stokes,  thirteen-moment.  Wang  Chang,  Burnett^and  Zoller 


solutions  all  tend  to  agree  with  one  anotlier.  An  attempt  has  recently 
been  made  (Haviland  )  to  solve  the  Boltzmann  equation  by  a  Monte 
Carlo  method  in  which  the  motion  of  a  statistically  large  number  of 
molecules  is  studied  individually.  The  procedure,  which  requires  a 
large  scale  computer  and  which  according  to  it.s  author  should  only  be 


applied  as  a  last  resort,  gave  results  for  shock-wave  structure  which 


seem  to  agree,  most  nearly  with  the  theoretical  calculation.s  according 
to  reference  39. 


In  appraising  the  Navicr-Stokes  equations  it  should  be  noted 
that  in  contra, St  to  solutions  along  stnmiy  kinetic  theory  linc.s,  the 
Navier-Stokes  equations  have  considerable  flexibility.  For  exampl.,, 
intermolccular  forces  can  be  easily  taken  into  account  by  permitting 
It  and  k  to  vary  m  a  prescribed  manner.  Moreover,  ,he  ratio  of  specific 
Heats  V  can  be  prescribed.  Finally,  it  is  even  possible,  in  a  compara- 
livMv  .-simple  manner,  to  take  into  account  relaivation  effects  for  a 


diatomic  gas,  In  appl^dng  the  Navier -Stokes  equations  (l}-(3)  for  a  dia¬ 
tomic  gas,  it  was  originally  found  (c.g,,  references  20,  22,  24)  that:  even 
for  fairly  weak  shock  waves,  the  observed  sliock  thicknesses  exceeded 
the  calculated  ones,  whereas  the  agreement  was  definitely  better  for 
monatomic  gases,  'i'liis  appears  to  be  due  primarily  to  the  rotational 
degree  of  freedom  of  a  diatomic  molecule  and  to  the  fact  that  the  attain¬ 
ment  of  thermal  equilibrium  for  the  rotational  degree  of  freedom  tends 
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to  lag  behind  the  equilibriirm  for  the  translational  degree  of  freedom 

42 

Gilbarg  and  Paolucci  have  pointed  out  that  a  thicker  shock  wave 
results  from  the  Navier -Stokes  equations  if  the  latter  are  modified, 
for  a  diatomic  gas,  by  introducing  a  bulk  viscosity  k  defined  by 
X  =  (2/3)(i+\,  where  \  is  the  compression  viscosity  (cf.  reference  2j  . 
Actually,  t  may  be  written  as  (cf.,  reference  2} 

T  =  -( 2|a  +  X)  u  '  (  35  a) 

InEqs.  (1)  -  (3)  and  in  Eqs  (34),  Stokes*  relation. 2ij  1  3X  =  0,  or 
K  -  0,  has  been  used,  if,  for  n  diatumic  gas,  this  relation  is  dropped, 
then  more  generally  X  =  x  -(2/'3)ii,  and 

-[(4/3)  t  x]u'  (35b) 

2 1 

For  air,  Gilbarg  and  Paolucci  and  also  Sherman”'  assumed  x=(2/3)u 
(corresponding  to  X  =  0)  ,  It  can  be  Seen  from  Eqs.  (33),  (34),  and  (35b) 
that  the  effect  of  thus  introducing  a  bulk  viscosity  is  simply  to  replace 
u  by  [U  T  (3/4)x],  i.c.,  the  effect  is  the  same  as  an  increase  in  the  shear 
viscosity  coefficient  p.  According,  for  example,  to  Wang,  Chang  and 
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Uhlenbeck  .  if  the  relaxation  times  are  sufficiently  small,,  then  (h/ij) 

can  be  related  to  the  relaxation  time,  Sherman^^  has  found  that  by  this 
means  the  Navier-Stokes  equations  lead  to  good  agreement  with  experi¬ 
mental  observations  for  cl..  a1  oruic  gases  for  kfach  numbei'S  up  to  2  0. 

A  rather  clear  comparison  of  experimental  data  with  results  of 
various  theories  can  be  found  in  references  21  and  27  \vhere,  briefly, 
it  is  concluded  that  for  Mach  numbers  up  to  about  2  0,  the  Navier- 
Stokes  equations  (modified  by  H  fordiaio  '.ic  gases)  1  ea.d  to  satisfactory 
agreement  with  experimental  obseryations,  while  none  of  the  other 
theories  leads  lo  any  better  agreement  in  this  Macli  number  range, 
with  several  of  them  (such  as  the  thirteen-moment  approximation)  leading 
to  poorer  agreement.  For  higher  Machuiumbcrs ,  it  would  at  present, 
perhaps,  appear  somewhat  premature  to  draw  any  definite  conclusions 
The  greatest  emmediate  need  in  this  area  is  additional  reliable  experi¬ 
mental  data  in  the  higher  Mach  number  range  and  a  reliable  mathematical 
estimate  of  tiie  magnitude  of  the  error  in  (he  more  promising  approxi¬ 
mate  solutions  (such  as  the  types  in  references  36  and  39)  of  tlie 
Bolt/.mann  equation  for  this  range  Of  course  still  more  accuraic 
.solutions  of  the  Boltxmann  equation  would  also  be  desirable 
Heat  Conduction Wiihoui  Viscosity  (Pr=:0) 

For  the  case  ot  u  =  0  (k/O).  Eqs  (lO)  and  (12)  yield  the  following 
parabolic  relauon  between  T  and  u 

T/T  =  {u/u^)[l  +•  vM/'(!  -  fu/u^))]  (36) 


Moreover,  Eq.  (1!3)  in  c onjunction  with  Eq.  (26)  reduces  to 


/(u  -  fi)  =  [mR(Y+  I ) /4k  (y -1)1  f  u  -  u  }(u-u  )  (3?) 


where 


8  =  u  (1  +  yM,^)/(2yM  2) 

From  Eqs.  (36)  and  (37) 

kT'  =  (m(Y+  1)/[2(y-  l)])(vk  -u)(u-u  )  (38) 

1  a 

Eq.  (37)  shows  that  as  u  varies  continuously  from  u  to  u  ,  cannot 

1  2 

vanish  except  at  the  end.s  (if  u/B  there),  Eq,  (38)  then  shows  that  T' 
will  be  zero  only  at  the  front  and  rear  ends  of  the  shock,  but  will  not 
be  zeJo  anywhere  inside.  In  particular,  T  ^  0  everywhere  within  the 

shock,  i.e.,  T  will  inc  I'ease  monotonically .  However,  according  to 
Eq.  (36),  dT/du  =  0  at  u  =  H.  Consequently  if  the  shock  is  strong  enough 
So  that  6^u^,  then  dh/du  would  vanish  for  some  value  of  u  inside  of  the 
shock  structure,  indicating  that  T  would  have  a  maximum  there  (cf. 

Fig,  3),  ill  contradiction  to  the  previous  conclusion.  Thi.s  means  that 
a  linnsition  from  state  i  to  state  2  witli  u  changing  continuously  across 
the  entire  shock,  would  be  impossible  in  such  a  case.  (T'hi.s  can  also 
bo  .seen  by  noting  that  according  to  Eq,  (37),  dx/du  =  0  at  u  =-.B;  if  thi.s 
occurred  within  the  shock,  then  there  would  result  two  different  values 
of  u  for  the  same  value  of  x).  Thus,  it  i.s  seen  that  a  completely  con¬ 
tinuous  shock  structure,  in  all  the  variables,  v.ull  occur  if  and  only  if 
the  shock  is  sufficiently  Weak. so  that  This  condition  (using 


Eqs,  (21))  can  be  written  in  any  of  the  i.ollowing  alternative  forms 


u 

2 

U 


1 

2yM,^ 


3y-1 

Y(3-y) 


u 
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>— 

3y-1 


(39) 


If  Y  =  I.4,  for  example,  condition  (39)  becomes  <  1.196.  A  result 
of  this  type  was  originally  found  by  Rayleigh^  and  it  indicates  an  inter¬ 
esting  difference  between  the  action  of  viscosity  and  heat  conduction. 

As  alieady  seen,  viscosity  alone,  without  heat  conduction,  is  capable  of 
supporting  a  steady  continuous  shock  structure  regardless  of  the  sliock 
sLiength  (cf,  Eq  (2?)),  on  the  oilier  hand,  heat  conduction  alone,  without 
viscosity,  can  siipport  a  steady  conti!iuou.s  sliork  structure  only  if  the 
shock  is  not  loo  strong. 

When  condition  (39)  is  satisfied  the  solution  of  Eq.  (37)  for 
eonstanl  k  is*"^ 
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(9  0) 


in  Lq  (40)  the  origin  lor  x  ha.s  been  chosen  so  that  u  ~  (u  -t-  u.^)/2  there. 
It  can  be  verihed,  by  solving  for  u'  in  terms  of  u  according  to  Eq  (37) 
and  di  llereutiaung  ihal  the  point  where  u  ''=0  will  be  different  from  the 
point  i;  L  t  (In  fact  when  M  ;s  slightly  ijclow 


'  necti  e,Kj)l;(. 

'  '  .>v  r  .)  ■ 

'•  '  re  A  and  B  a; 

constant 


i.'vi’d  that  in  ati  lour  cases  treated  here  for  which  x  vs, 
iiiy  fourd  [Pr=3/4,  0,  o;i  with  constant  u  and  k;  weak 
solutions  arc  ol  the  Jorm  x-A  Iog(u; -u)-i3  log(u-U;))-t  C , 
t  , c  e !  t <  1 ! n  CO n s t a n i  s ,  a n d  C,  i $  a n  is  r b 1 1  r b  ry  i n t. e g r a i  i o n 
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(3y-  1)/[v{3-  y)]j  the  inflection  point  on  the  velocity  curve  will  be  far 
downstream,  nearu=Ug). 

If  u  /u,  <(y+  1)/(3Y"'  1)j  it  has  already  been  seen  that  a  continuous 

solution  for  the  entire  shock  structure  will  not  exist.  To  see  the  nature 

of  the  solution  in  this  case,  it  is  recalled  that  since  8>u_^  now,  and 

dT/du  =  0  at  u  =  8,  the  curve  T  vs.  u  will  have  a  maximum  between  u 

and  u  =  (see  Fig.  3).  This  means  that  in  this  case  the  temperature 

will  reach  its  final  value  T  =T,  at  a  value  of  u,  u=u^  (say),  between 

and  u^.  An  isothermal  discontinuity  then  occurs  at  this  point,  the 

temperature  remaining  at  T  but  the  velocity  (and  also  other 

variables)  jumping  discontinuously  to  its  final  value  u^  (c£.  reference 

14).  The  actual  T  vs.  u  curve  is  thus  the  curve  lc2  in  Fig.  3,  con- 

1  2 

si  sting  of  the  arc  Ic  and  the  horizontal  line  segment  c2.  Gilbarg  has 
proven  malhcmatidally  that  as  p  approaches  zero,  the  solution  of  Eqs, 

(13)  and  (14)  approache.s  this  type  of  discontinuous  solution  Tlic  location 
of  the  point  "c"  in  Fig.  3  can  be  found  by  .setting  T-T^  in  Eq.  (36)  anti 
solving  for  u.  One  root  must  be  u  and  the  oilier  root  will  be  u  =  u^ 
w  h  e  r  e 

u^/u^  =  [2yM/  -  (v-I)]/[v(Yf  1)M^"]  (41) 

The  velocity  u  will  (iet  rea.se  continuouBly  with  from  at  x  =  -  ix>  until 
it  reache.s  the  value  ,  at  which  point  it  then  drops  cii  scontinuou.sly  to 
u^  and  JU'inain.s  ai  'h.al  I'alue.  Ciirvt-.s  showing  T  and  u  v.s,  x  in  both  the 
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This  IS  not  t  rue  of  the  density  or  p.rcssurc,  which  wi.ll  vary  mono- 
toirirally  with  u  (cf,  Eqs.  (6)  and  (10)  with  l  =  Q), 


The 


continuous  and  discontinuous  cases  are  given  in  Figs,  4  and  5. 
special  case  of  =  (3 y- l)/f  Y(3- y)]  is  also  included  there. 

IV  -  DISTRIBUTION  OF  ENTROPY 

It  has  already^  been  seen  that  the  velocity,  specific  volunie,  pressur 
and  temperature  all  vary  monotonically  across  the  entire  shock  width, 
the  first  two  decreasing  and  the  latter  two  increasing  throughout  (with  the 
indicated  modifications  for  Pr^rO).  The  entropy,  however,  although  it 
will  be  greater  at  the  end  of  the  shock  (x  =  +  oo)  than  ahead  of  it  (x^-o)}, 
will  be  iound  not  to  increase  throughout  the  shock  structure,  but  to 
increase  from  x--oo  until  it  reaches  a  maximum  somewherewithin  the 
shock,  and  then  to  decrease  until  its  final  value  at  x  =  +oo.  Thi.s 
anomalous  behavior  of  the  enli'opy  will  now  be  examined  in  detail,  based 
on  the  Navier-Stokes  equations,  first  for  general  flows  [i.e.,  general 
Prandtl  number),  and  then  for  tlic  specific  flows  in  Section  Hi  foi' 
which  exact  closed-form  solutions  have  been  found.  For  convenience 
the  results  obtained  here  will  be  pul  in  the  form  of  numbered  theorems 
and  corollaries  and  these,  can  then  serve  as  a  convenient  summary  of 
the.  results  developed  here, 

Gene ral  Flows  (Pr  general) 

III  t:.i'  tir.-'’  pl.ji  (•  It  hi’  th.it  tile  i.i’iM  Luum.i  litut  iif 


1‘ '  c  ‘1  iil  fiicicru  «•  sK  I  1 1  )•  1  •;  ti  .  x  :n  such  .i  i.m*’  is  not  auilc 

a.  I  u  .-.it  .  -itu  I-  It  .i-.fii,  .iti’S  li-.it  1.  m;11  v.i  ry  <  out  i  r;- 1 1  ■ , I V  .iiijl  it  ha.s  the 
\;.li;r  ;  {•Ai’,1  rr  r  -  I'llmile  nd  th>’n  rirnp^  abruptly  tous. 

.\ e  1  u.ii  I’.- ,  i.-T  v.iricr,  ^  I'ru .  i;i.o  .  s  1  \  only  i"i’l  l‘,,iM  the  v.il'.i'  ^  which  will 
Oil  -i'  uiiHin  ,im  ;n,'  po;;::  w},i  vr  T  '  -  d  (  i  F 1 1;  A  t;fr.),  .md  then  drops 
'u*  I'l.i  lett  ‘'.indi  r!  fit’ ,  .  .w.  i .  be  1 .  iii:  e  -t '  (cf.  Fig.  Shore). 


/  ’4 


•  i  - 


the  gas  plays  a  key  role  in  the  non- monotonic  behavior  ot  the  entropy. 

This  can  be  readily  seen  by  noting  that  for  a  gas  witliout  heat  conductioTi, 
i,  e.  ,  for  a  gas  with  k  =  0,  the  energy  equation  would  imply  S  0  within  « 
the  entire  shock  structure;  hence  in  this  case  the  entropy  would  indeed 
increase  th roughout  the  shock.  Physically  this  is  simply  due  to  the  fact 
that  heat  keeps  being  added  to  the  gas  due  to  the  viscous  dissipation, 
while  there  is  no  heat  being  conducted  away,  Thus, 

Theorem  I,  In  the  absence  of  heat  conduction,  but  in  the  presence  of 
viscosity,  the  entropy  would  increase  rnonolonically  ihrotighoul  the 
entire  shock  structure, 

In  the  presence  of  heat  conduction  (k  >0),  however,,  it  is  actually 
not  difficult  to  see  that,  in  genei'al,  ihe.  entropy  can  no  longer  vary 
rnonolonically.  Tliis  follows  by  considering  condition.s  dost'  to  both 
ends  (x=±oo)  ol  the  sltock  wave  Because  of  the  strongly  vmiform  cbncli- 
iions  near  cacli  end,  tlie  derivative  will  be  very  small  in  tlmse  vioimtics, 
and  consequently  Uu '■"  will  be  much  smaller  than  (kT')'  tiicre.  llenCc, 
according  to  Eq  (9),  S'  wnll  have  the  same  sign  as  (kT')',  or  T"  {.since 
k  ^T  may  also  be  neglected  here  and  k>0)  in  these  neighborhoods.  Since 
T  vs,  X  will  in  general  follow  an  S  curve  (c(.  ,  e.g,  ,  Fig,  Z),  T">0 
near  x  =  -cxi  and  T^'<0  near  x  r  f  co-  Con.scquent.ly  S  must  be  inc  reasing 
near  the  upstream  end  ol  the  shock  structure,  bul  must  be  decreasing 
non  r  the  down  St  ream  end 

"  It  is  intere.sling  to  note  th.it  despite  the  as-ipareni  .simplicity  of  this  argu¬ 
ment,  the  very  first  reaction  oi  most  investigators  in  this  or  related  fields 
on  learning  of  this  fact  has  been  that  of  surpri.se. 


The  above  argument, which  gives  soarie  insight  into  tire  cause 
(namely,  the  role  of  heat  conduction)  of  the  anomalous  behavior  of  the 
entropy, can  be  made  mathematically  more  formal  by  considering  solu¬ 
tions  in  the  (u,T)  plane,  as  in  reference  16.  Multiplying  Eq.  (9)  by  dx, 
mTdS=d(kT  )+  [4/3)lJuMu.  Hence,  from  Eqs.  (13)  and  (14) 


niTdS  =  dL  +  (4/3)  Mdxi 


(42) 


At  the  end  points  I  and  2  of  the  shock  wave,  L  =  0  and  M  =  0  (cf.  footnote 
3).  Hence,  at  each  end  point  rnTdS^dL.  However,  L>0  everywhere 
within  the  shock  structure  (references  11  and  12).  dL>0,  and  hence 

dS>0  at  the  front  end,  but  dL<0  and  hence  dS  <  0  at  the  rear  end.  Thus, 
Theorem  2  follows. 

Theorem  2.  In  the  presence  of  heat  conduction,  the  entropy  will  not 
inciease  monotonical ly  throughout  the  shock  structure,  but  will  increase 
rmar  the  upstream  edge  of  tiie  wave,  and  if  it  varies  continuously,  will 
decrease  near  the  downstream  edge. 

Corollary.  If  k>0  and  S'  is  everywhere  continuous,  then  at  some 

point  within  th.e  shock  wave. 


In  (.  onncctioti  with  Iheorem  2  ar.d  its  Corollary,  it  should  be 
recalled  that  S  and  will  in  fact  vary  continuously  in  alt  ta,se3  e.xcept 
lliai  0.1  a  suJficiently  strong  siiock  v.nve  wjthou!  visco.siiv.  It  should, 
of  comse,  be  immcdtaleiv  noted  here  that  Theorem  2  does  not  violate 


the  second  law  of  thermodynamics,. 


since  Only  the  entire  shock  strnriure 


with  its  uniform,  conditions  a«  x  =  ± 
hence  it  Ls  only  required  that  S„  >S 


00  con.stitutes  a  closed  .system,  and 
,  .  Further  remarks  on  the  physics 


of  Tlieorem  2  will  be  n  ade  subsequeiitly , 

From  Eq,  (9)  it  is  s een  that  S 0  wh ere  (kT ^  ,  or 
where  =  - (4/3  )|iu 'T Hence,  ifk>0,  jj  >  0  and  k '>  0  throughout 
the  shock  (as  is  ordinarily  the  case),  then  S  '=0  at  a  point  within  the  shock 
where  T^^<0.  Keeping  in  mind  the  general  S-shaped  nature  of  the  tempera 
ture  distribution  (cf.  ,  e.  g.  ,  Fig.  2),  T"<0  downstream  of  the  point  where 
T"=0.  Thus,  Theorem  3  follows, 

Theorem  3,  If  both  heat  conduction  and  viscosity  are  present  and  the 
heat  conductivity  coefficient  is  either  constant  or  increasing  throughout 
the  shock,  then  the  entropy  will  have  its  maximum  value  at.  a  point  within 
the  shock  structure  which  is  downstream  of  the  point  of  inflection  of  th ^ 
tcmperaturc-vs.  -distance  curve. 

If  Eq,  (2)  is  multiplied  by  (dx/dp)  and  use  made  of  Eq,  (6),  thc'n 
it  is  found  that  (dp/d  p)-u^ -(4/3)(uu ')  ^{dx/dp)=  0.  Tliis  equation  iiolds 
everywliere  inside  of  the  shock  where  do/dx/O  and  hence  at  a  point  inside 
of  the  shock  where  S  '  =  0,  Since  at  such  a  point  (dp/dp}  =  a'^,  Theorem  4 
f  ol  1  o\v  s  , 

Tlieorem  4.  At  anv  point  inside  of  the  .shock  wave  where  S'-O, 
u^^a--(4/3)(dx/dD)(!Ju')^ 

This  theorem  yields  the  foil owi ng  impo riant  corollary, 

Corollary.  1  f  at  a  point  inside  of  tlie  sliock  wave  w-here  S  '  -  0  it  is  true 
that  (lau')'- 0,  then  in  addition  u=a  there,  i,e.  ,  the  speed  must  be 
locally  sonic  there  Conversely,  if  at  a  point  inside  of  the  shock  where 


S'=Q  it  is  true  that  the  velocity  is  locally  soniC;  then  also  (pu^)^=Q  there. 

It  should  be  noted  that  in  the  cases  of  a  constant  non-zero  viscosity 
coefficient  [J,  the  condition  (uu^)^=0  means  an  inflection  point  in  the  vclocity- 
vs.  -distance  curve.  Details  of  the  entropy  distribution  for  various  special 
cases  will  now  be  analyzed. 

Weak  Shock  Waves 


For  this  case  it  will  be  recalled  li  and  k  may  be  considered  constant, 
while  is  second-order  small,  Consequently  the  term  may  be 
neglected  in  Eq.  (9).  Moreover,  since  u''  is  third-order  small,  so  is  T" 
and  hence  one  may  replace  T  by  on  the  left  side  of  Eq.  (9),  solve  for 
S  '  and  integrate  to  obtain  S-S^  =  [k/(mT^)]T  .  From  Eq.  (31)  it  then 
follows  that 


S-S^  = - (T^-Tjtu  -u  )  sech^  x  (43) 

*  4mDT^  -  I 


A  result  analogous  to  Eq.  (43)  has  been  obtained  by  a  somewhat  diffcrc'nt 
procedure  by  Landau  and  Lifshitz^"^,  According  to  Ea.  (43)  S-S^  varie.s 
symmetrically  about  x  =  0  and  has  its  maximum  value  there  (see  Fig.  6). 
Ihe  point  x  =  0  is  also  where  u,  p,  o  and  T  each  have  a  value  equal  to 


"In  this  connection  it  must  be  remarked  that  Golit.syn  and  Staniukovich 
have  .stated  (for  constant  u)  tliai  in  general  u"=0  wlmre  S''=0,  and  hence 
have  concluded  also  that  in  gencrar’s'  =  0  where  u=a.  Although,  a.s  will  be 
seen  for  corivStant  y,  a'^=0  at  =  0  for  weak  shock  waves  and  for  any  .shock 
waves  when  Pr=3/4,  Golitsyn  and  Staniukovich  do  not  give  any  proof  nor 
any  reference  lor  the  truth  of  this  statement  in  general-  (They  do  give  a 
reference  for  the  proof  of  this  statement  for  weak  shocks.)  This  state¬ 
ment  of  Golitsyn  and  Staniukovich  must  at  present,  therefore,  remain 
open  to  question  for  the  general  case  when  u/O.  When  U=0,  it  will  be 
seen  that  the  statement  that  u'^^O  where  S'^O,  is  false. 


the  arithinetic  mean  of  their  values  in  front  of  and  behind  the  shock  (cf. 
Eqs,  (30)  and  (31)).  Moreover,  at  this  point,  according  to  Eqs.  (30)  and 
(31),  u",  p",  p^',  T"  all  vanish.  Hence,  from  Theorem  4,  the  velocity 
is  locally  sonic  there.  These  results  can  be  summarized  by  the  following 
theorem. 

Theorem  5,  For  extremely  weak  shock  waves  (vath  k>  0),  the  entropy 
has  its  maximum  at  a  point  about  which  it  varies  symmetrically.  At 
this  point,  the  velocity  is  locally  sonic.  Moreover,  the  velocity,  dejisity, 
pressure  and  temperature  curves  vs,  distance  all  have  their  points  of 
inflection  there,  while  each  of  these  variables  has  a  vaUic  equal  lo  the 
arithmetic  mean  of  its  values  ahead  of  and  behind  the  shock. 

[t  will  be  seen  that  for  the  other  cases  to  be  analyzed  !u:re,  some 
of  these  statements  continue  to  hold,  while  others  do  not..  It  may  be 
noted  (hal  setting  x  =.  +  co  in  Eq.  (43)  it  is  found  that  So-S-,  .  This  is  a 
reflection  of  the  fact  that  for  extremely  weak  shocks,  is 

actually  only  third- ordc r  Small  (cf,  reference  14  or  44);  the  maximum 
value  of  S-Si  within  the  shock  (x~0)  is,  however,  one,  order  of  magni¬ 
tude  larger,  namely,  second- order  small  This  limiting  case  of  weak 
shocks,  therefore,  reveals  expecially  strongly  the  non-monotonic  character 
of  the  entropy  variat,ion  through  the  shock, 

Prandtl  number  of  3^4 

From  ,Eci,  (331,  v.hich  Isolds  for  lliis  case.  T'  ~  uu'/c  ,  Substi- 
luting  for  T.'  into  Fq.  (9),  and  also  .substituting  k:;{4/3)uc^  there,  it  is 

cTS'  ::  -(4/3)(au*)  '  (44! 

‘j  4 


found  that 


Integrating  Eq.  (44)  over  the  entire  shock  and  noting  that  u^=0  at  each 
end,  one  obtains 


Z  2 

[  pTdS  =  0  ,  or  f  pdS  =  0  (45) 

a  J, 


Since  p  is  everywhere  positive,  the  integral  in  Eq.  (45)  can  vanish  only 
if  S'"  is  positive  in  some  places  and  negative  in  others.  Hence,  S  cannot 
increase  monotonically  across  the,  entire  shock.  This  result  has,  of 
course,  already  been  proven  for  more  general  cases  (general  Prjifco), 
but  the  proof  given  here  is  interesting  because  of  the  rather  ijileresting 
result  given  by  Eq.  (45)  and  because  this  is  the  reasoning  first  used  in 
1944  by  Eoy^,  who  may  have  been  the  first  to  discover  the  non-monotonic 
behavior  of  tlie  entropy.'-®  Moreover,  equations  analogous  to  (44)  and  (45), 
but  with  heat  conduction  replacing  the  role  of  the  viscous  shear  stresse.s 
here,  will  subsequently  be  found  to  hold  for  a  fluid  with  heat  conductiou  but 
fi-oi't  viscosity  (P.r:;0). 

From  Eq.  (44)  it  follows  that.  S '  =  0  where  (mu')'  =  0.  From  Eq.(Z.i) 
it  is  recalled  that  this  is  the  point  where  u=/^/u,Ug  ,  Moreover,  since 
Liu  '  denotes  the  shear  stress,  the  point  where  S  is  a  maximum  is  in  this 
case  also  the  point  where  the  shear  stress  is  a(  no  ga  tiv  e)  maxi  mum , 

Finally,  from  Theorem  4,  it  also  follows  that  u  =  a  there.  These  r<>su!ts 
all  hold  for  variable  a.s  well  as  constant  |a  and  can  be  summarized  by  uic 


1“ 

In  .hoy's  paper  the  factor  "c”  was 
Eq.  (45)  it  was  inferred  there  that. 


omitted  in  Eq. 
2 


TdS  1.  0. 


(4,4),  and  in. stead  of 


I 
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following  theorem. 

Theorem  6,  For  a  constant  Prandtl  nruTiber  of  3/4,  but  otherV'hse  variable 
viscosity  aoid  heat  conductivity  coefficients,  the  entropy  will  have  its 
maximum  value  within  the  shock  W'ave  at  the  point  where  the  shear  stress 
is  a  maximum.  At  this  point  the  velocity  is  locally  sonic  and  has  a  value 
equal  to  the  geometric  mean  of  the  values  of  the  velocity  ahead  of  and 
behind  the  shock 

Corollary  1,  If  Pr=  3 /4  and  the  viscosity  (and  therefore  the  heat  con- 

cUictivily)  coefficient  is  constant,  then  the  entropy  has  its  maximum  vahic 

at  the  inflection  point  of  the.  vclocity-vs,  -dislajice  curve. 

If  u  is  variable,  the  S=S  where  i.i  u'  -1  i',u'  TO.,  nonce,  sbreq 

max  ' 

u'  <0,  it  follows  that  if  [a  is  an  increasing  function  of  'i'(as  is  ordinarily 

the  case),  so  that  Li'  >0,  then  S  occurs  where  u''  >0,  and  therefore 

max 

downstream  of  the  inflection  point,  on  the  velocily-di st ance  cvirvc.  7'hus, 
Corollary  2.  If  Pr'.;3/4,  and  u(and  k)  varies  so  that  it  increases 
throughout  the  shock,  then  the  entropy  w-il!  have  its  m;ixi,mum  value  at  a 
point  downstream  of  the  inflection  point  on  the  vdocity-vs,  -distance 
curve , 

In  actual  cases  (cf  Fig  7)  the  point  of  maximum  >S  may  still  be 
close  to  the  inflection  point  on  tiio  velocity-distance  curve.  Typical  curves 


showing  ihiC  cm  ropy  distributloii  are  giv<.’n  hi  Fig,  7., 
can  be  .found  in  roferemces  9,  7,  and  15. 

Viscosily  VvCthout  Heat  Conduction  (Pr=a;)) 

It  has  already  beeii  poijtted  out  that  this  is  the 
the  LT7n.rQpy  v\i!!  increav^ic  monotonically  throughout  i!i 


V  d  d  i  t  X  o  n  a  1  c  u  r  v  ^  *  s 


oric  case  in  wlxicii 


shock  This  is 


iiiciSt  r,-ih:’f.i  xn  Fig,  8.  It  may  be  added  that  here  Is-a  case  tn  v/pich  the 
■  ■  'f  •  ■  also  the  .  ui  have  alrnosi,  cxactlv  as  for 


th*'  case  Pi-3/4  (compare  Eqs.  (27)  and  (23)),  wMle  the  enti’opy  has  a 
qualitatively  different  behavior. 

Heat  Conduction  Without  Viscosity  (Pi-Q) 

The  case  of  a  completely  continuous  shock  structure 
(M^  <[3y-1)/[y(3-y)])  will  be  discussed  first.  Putting  M  =  0  into  Eq .  (9), 

it  is  seen  that  in  this  case 


mTS'=  {kT')  ' 


(46) 


From  this  it  follows  that 

2 

•» 

TdS  =  0 


(4-7) 


r^qs.  (46)  and  (47)  are  analogous  to  Eqs.  (44)  and  (45)  for  Pr=3/4. 

Eq,  (47)  indicates  that  vvitliout  viscosity,  the  net  heat  added  to  the  fluid 
as  u  crosses  the  shock  is  zero,  i.  e.  ,  without  viscous  dissipation,  there 


will  be  as  much  heat  conducted  away  from  the  fluid  as  into  it. 

From  Eq.  (46)  it  follows  that  S  0  where  (kT  ']  '=0.  Hence, 


^"^max  where  the  heat  flow  kT  '  is  a  maximum.  More.over, 

differentiating  the  right  side  of  Eq.  (38),  it  is  readily  seen  that  (kT')'=0 
vhcie  u~(u^t  u^)/2.  It  has  been  previously  noted  that  at  tliis  point, 
(except  in  the  limit  of  weak  shocks).  Since  the  pressure  p  in  this  ca.se 


vanes  linearly  with  u  (Eq,  (10)  with  u-O),  it  follow.^  also  that  ps(p  fp^)/2 
ai  this  point.,  (This,  however,  will  nm  be  true  of  the  density  or  tempem- 
tuic),  Finally,  .since  u  =  0  here,  Theorem  4  show.s  that  u-n  a*  thi.s  point 
These  results  all  hold  for  variable,  as  well  as  constant,  k  and  ran  be 


summarized  by  the  following  ihcorem 


Theorem  7.  If  heat  conduction  only  is  present,  without  viscosity,  and 

tVi  P  S^Qf"  V  \r  ic  cnf-Tir’nriv^fltf  e.r>  fK-ni-  ^  4-1,  q,-.  fl,  , 

v^v<..<Aiv^xv,nvi.y  i.jic3.'w  i  Vt .  ^  \~^Y~-^^/[.  yw  Y/J)  11 1  c  u  IJ  m 

entropy  will  vary  continuously.  It  will  reach  its  maximum  value  where 
the  heat  flow  is  a  maximum.  At  this  point,  the  velocity  is  locally  sonic 
and  both  the  velocity  and  the  pressure  are  equal  to  lire  arithmetic  mean 
of  their  values  ahead  of  and  behind  the  shock. 


F rom  the  fact  that  S  ^  =  0  where  (kT  ')  '=  0,  the  corollary  below 

follows. 

Corollary  3.  Under  the  same  conditions  as  in  Tlieorem  7,  the  point  at 
which  the  entropy  reaches  its  maximum  will  be  at  the  point  of  infl  ecti on 
of  the  tempo ratii re- vs ■  -di.stance  curve  if  the  heat-conductivity  coefficient 
k  is  a  consUnl,  and  will  be  downstream  of  this  inflection  point  if  k 
increases  through  the  shock. 


This  corollary  may  be  compared  with  Theorem  3. 

If  >  (3  y- 1 )/[  y(3  -  v)],  It  will  be  recalled  tliat  the  shock  structure 
is  continuous  until  T=T  ,  whereh  =  u  .  At  this  point  an  isotliermal  dis- 

2  C 

continuilv  occurs.  The  entropy  in  such  a  case  will  vary  contintiou.sly 

froin  point  "!*'  to  tlii.s  poi.nt  "c”,  and  will  th.cn  change  abruptly  (at  point 

b  )  frorn  its  value  there  to  its- final  value  S..^ ,  Tv.-'c)  (luesltons  of 

particular  interc.st  liert'  ai'C'  (a)  flow  does  S_  compare  with  S  .  i,  c.  , 

c  ?  ■ 

wlien  the  entropy  jumps  from  lo'S^,  does  it  jump  to  a  higher  or  to  a 

lower  value?  (b)  During  tlu;  cqniinuous  variation  of  llu;  entropy  from 

^5  ‘  docs  t'ne  entroj)y  iticrca,sc  monotOnicaily -to 

S  ,  or  does  it  at  first  increase  to  some  value,  S  .  wfierr  it  k;,,-.-  a 
^  '  m,ax'  “  ’  ' 

local  maxi  tnum  witn  zero  slope  !S^~0)  and  then  dee  re?,  .ue  to  -S  ^ 


Question  (a)  is  readily  answered.  Since  T^  =  T^  and  u^>u^.  it 
follows  from  Eq,  (8)  that  S^>S^.  Hence,  in  the  isothermal  discontinuity 
here,  the  entropy  jumps  to  a  lower  value.  This  is  thus  a  (comparatively 
rare)  example  of  a  physical  system  within  which  a  discontinuity  to  a 
lower  entropy  occurs.  In  regard  to  question  (b),  S  will  pass  through  a 
maximum  with  zero  slope  before  the  discontinuity  if  (u^ +u  )/2>u  . 
Substituting  for  u^  and  u^  in  terms  of  M  according  to  Eqs,  (21)  and  (41), 
this  condition  becomes 

<(2y-1)/[v(2-y}] 

For  y=J-4,  for  example,  this  condition  is  M^<  1.460.  These  results 
for  the  discontinuous  case  can  be  .summarized  by  the  following  theorem. 
Theorem  8.  If  heat  conduction  only  is  present,  without  viscosity,  and  the 
shock  is  sufficiently  strong  so  that  M ^ ^  >(3 Y~  1)/[y(3- y)1,  then,  the  entropy 
will  vary  continuously  within  the  shock  strvicture  until  the  temperature 
first  attains  its  final  value  T  .  At  (his  point  an  isothermal  discontin- 

- -  -I . . . . . ———I  — . .  45  -  '  “ 

uity  occurs  to  the  final  state  behind  the  shock,  with  tlie  cnlro];y  drop]?ing 
abruptly  from  a  higher  value,  here  to  a  lower  value.  If,  in  pa)'iiciiiar , 
(3y-1)/[y[3-v)K  M ^ ^  <(2y- 1)/[y(2-y)],  then  the  entropy  will  fi rst  inc reasc 
continuously  until  it  reaches  a  maximum  value  with  zero  slope  (v.di o re 
u=(u.  1-u  )/2=a),  decreases  continuou.s]y  until  T-T  ,  and  then  drops 
a bruptly  to  its  still  lowe r  fi nal  value  S ^  Finally,  if  M. ^>(2y-  1  )/I v( <1- y)]. 
then  the  enironv  will  increase  continuouslv  until  T  =  T  ,  and  then  vhll 

. - . . .  . . . . .  I  '■  "  ^  — — I— ■  '■ 

drop  abruptly  to  its  final  value. 

Theorem  S  (a.s  svell  a.s  Theorem  7)  hold.s  for  variab'e  a.s  well  as 


constant  k,  Fig.  9  shows  various  types  of  entropy  distribution  in  tlie  case 
of  heat  conduction  without  viscosity, 

Physical  Remarks 

The  mathematical  results  developed  here  on  the  entropy  distribu¬ 
tion  warrant  at  least  some  physical  comment.  From  a  macroscopic, 
thermodynamic  point  of  view,  the  physical  explanation  for  the  non¬ 
monotonic  distribution  of  the  entropy,  as  already  indicated,  is  relatively 
simple  and  is  based  on  the  heat  conductivity  of  tlie  gas.  The  licat  added 
per  unit  time  to  a  unit  volume  of  tlie  gas  at  any  point  is  (kT  ')  ^  Because 
of  the  nature  of  the  temperature  distribution  (e.  g,  ,  ]''ig,  Z)  wi(liin  Ihe 
shock  structure,  (kT ^  is  positive  in  the  front  part  of  the  shock 
structure,  and  negative  in  the  rear  part.  Consequently,  in  the  rear  part 
the  fluid  loses  heat  through  conduction  and  when  the  tiine  rate  of  this  loss 
exceeds  tlie  rate  of  gain  ol  heat  by  viscous  dissipation,  the  entropy  (as 
defined  liere,  viz.  dS  =  dQ/T)  will  dimmish.  This  is,  of  course,  the 
content  of  Theorems  1-3  and  their  proofs. 

The  foregoing  .simple  explanation  for  the  behavior  of  tlie  eniropy, 
however,  docs  not  leave  one  entirely  satisfied  physically,  since  if  is 
ba.scd  on  considering  entropy  change  to  be  simply  proporiional  lo  lieai 
added  (albeit  with  a  variable  factor  i/T)  and  is  not  concerned  w'ith  the. 
.special  si  gni  ficance  of  entropy  it.se!  f,  namely,  a.s  a  measure  of  the 
"disorder*'  of  a  sy.stcm.  (In  statistical  meclianics,  a  system  is  associated 
with  large  "disorder"  or  "randomness"  if  its  microscopic  jiroperti  e  s , 
e,  g.  ,  the  distribution  of  molecular  velocities  for  a  gas,  may  be  arranged 
in  very  many  ways,  all  consistent  with  the  same  macroscopic  propcrtic,? 
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(e,  g.  ,  the  temperature  and  pressure  for  a  gas)  of  the  system  {cf.  , 
e.  g,  ,  reference  45)).  The  results  developed  here  indicate  then  that,  in 
the  presence  of  heat  conduction,  the  state  of  the  gas  throughout  a  shock 
wave  varies  in  such  a  wav  that  although  the  state  of  disoidei  iS  gieatei 
behind  than  in  front  of  the  shock  (second  law  of  thermodynamics),  there 
is  an  even  greater  state  of  disorder  at  certain  points  within  the  shock 
St  rue  tur  e . 

The  chief  difficulty  in  giving  a  completely  satisfactory  statistical 
mechanical  interpretation  of  the  results  for  entropy  obtained  here  is  that 
in  a  shock  structure,  there  are  large  deviations  from  statistical  equi¬ 
librium,  i.  e,  ,  there  are  very  large  gradients  of  the  macroscopic  .state 
of  the  gas  within  the  shock,  In  such  a  ca.se  the  definition  of  entropy 
S  as  used  here  should  be  replaced  by  the  more  general  lunclion  H/r, 
where  H  is  Boltzmann's  H-function  in  tlie  kinetic  theory  o!  gases,  1'  or 
Maxwellian  distributions  of  velocity,  H/p  is  Indeed  .simply  proportional 
to  the  negative  of  S,  but  for  appreciable  deviations  from  equilibrium, 
no  siich  .simple  relation  between  S  and  H/p  exists.  (See,  c.  g-  ,  refei- 
encc.  11  lor  furtlicr  cleinils  and  derivations).  It  ir.ay  still  be  expee  1  ed , 
liowt'ver,  tha.,1  at  least  for  .chiH:  i,  vv.ives  noi  ido  strong,  the  entropy  S 
should  correspond  fairly  well  to  H/p.  It  is  noteworthy  tliat  in  all  of  the 

Another  meaning  of  large  deviation  trom  equilibrium  is  that  the  mole¬ 
cular  di  .St  ribution  of  velocities  is  far  from  MaxuvelHan  (cf.  ,  c.  g.  , 
reference  2'^''),  For  the  shock  structure,  the  distribution  of  molecular 
velocities  is  Maxwellian  in  the  uniform  conditiotis  ahead  of  and  behind 
the  shock,  but  will  not  be  ex.aclly  Maxwellian  within  tiie  .shock,  and  may 
be  expected  to  be  furthest  from  Ma.xwelHan  around  the'tentcr"  of  the 
w.ive 


•1] 


special  cases  analyzed  here  in  which  the  entropy  5  was  non-rnonotonic 


(k  /Q) ,  S  had  its  maximum  value,  i.e,  ,  the  ''disorder"  was  a  maximum) 

where  the  gradients  of  certain  physically  significant  quantities  (such  as 

the  shear  stress  or  tlie  heat  flux)  were  a  maximum  (and  perhaps  where  the 

deviations  from  a  Maxwellian,  distribution  might  be  larjrest).  It  may  also 

be  significant  that  inevery  case  analyzed  here;  the  velocity  at  this  point 

was  locally  sonic,  i.e.  ,  this  is  the  point  where  the  flow  changes  from 

supersonic  to  subsonic,  Further  interprcladoiis  along  these  lines  can  be 

made  only  from  reliably  accurate  solutions  of  the  Boltzmann  equation, 

and  calculations  of  H  therefrom  It  is  of  special  interest  to  note  tliat 
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Haviland  has  actually  obtained  distributions  of  H/p  throvigi:  a  shock  struc 
Cure  fur  M^=2.0and  3.0,  by  a  Monte  Carlo  procedure,  and  ha, s 
found  H/p  to  pass  through  a  minimum  (corresponding  to  a  maximum  of 
S)  in.side  of  the  shock  structure. 

V  -  STRUCTURE  OF  SHOCK  WAVES  INVOLVING  CHEMICAL 

reactions 

In  the  previous  sections  there  have  been  discu.ssed  the  structure 
of  shock  waves'  in  general  and  the  entropy  distribution  therein  in 
particular;  the  classical  assumptions  of  a  homogeneous  gas  with  constant 
coefficients- of  specific  heat  were  employed.  It  is  the  purpose  of  this 
section  to  review  the  analysis  applicable  to  the  structure  of  shock  wave.s 
which  involve  chcnucai  reactions  in  the  downstream  region.  Such 

^'The  analysis  in  this  section  is  a  review  and  generalization  of  ilie 
stuciic.s  of  Hi  rschfelcier ,  Curuas  and  coworkers  (cf.  ,  e.  reference  4  6) 
periaimng  to  detonation  waves  and  of'Gai  tatzes  and  Bloom‘‘~  pertaining  to 
sironR  .shock  waves  in  dia.m  'cgascs  vSee  ai.$o  the  rcCent  excellent 
re\dcw  of  deronati on  waves  by  Oppenh ei nt 'and  Stern"^^. 


reaction  arises  when  a  chemically  reactive  mixture  of  gases  upstream 
of  the  shock  are  raised  by  the  initial  compression  in  the  shock  to  a  state 
in  which  significant  chemical  reaction  takes  place.  The  steady  state 
conditions  across  the  waves  are  altered  by  the  chemical  reaction. 

This  type  of  wave,  when  the  mixture  upstream  is  a  mixture  of 
fuel  and  oxidizer,  is  called  a  detonation  wave;  it  has  been  the  object  of 
investigation  for  roughly  80  years.  Particular  attention  has  been  devoted 
to  detonation  waves  satisfying  the  Chapman- Jouguet  condition,  Uiat 
infinitely  far  downstream  the  flow  velocity  is  equal  to  the  local  sonic 
velocity,  It  has  been  long  recognized  that  this  condition  is  not  unique 
but  corresponds  to  the  physical  situation  usually  occurring  in  the  labor¬ 
atory  or  in  nature,  e.  g.  ,  in  the  case  of  mine  explosions.  Recent 

experimental  tecluiiques  employing  special  supcJ'sonic  wind  tunnels 
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(cl,  ,  e.  g.  ,  NichoUs  Gross  and  Chinnitz"  ,  and  Rhodes  et  al 

have  afforded  opportunities  to  study  detonation  waves  not  satisfying 

the  C  -  J  condition. 

It  is  al.so  possible  to  obtain  chemical  reactions  behind  shock 
waves  in  diatonkc  gases  such  as  air.  If  the  shock  is  sufficiently  strong, 
i.e.  ,  involves  temperatures  and  densities  in  its  downstream  region, 
sufficiently  high,  dissociation  of  the  gas  molecules  takes  place  and  the 
gas  composition  involves  both  atoms  and  molecules.  The  dissociation 
energy  of  the  molecule  is  absorbed  by' the  gas  mixture  altering  both 
velocity  and  temperature.  This  is  a  form  of  chemical  reaction  which 
has  been  under  intense  study  because  of  it.s  importance  in  hypcreonic 
aerodynamics.  Hovrever,  the  only  study  of  the  structure  of  .such  waves 


known  to  the  authors  is  that  of  Gaitatzes  and  Bloom  ,  In  general  the 
physical  situation  leading  to  the  establishment  of  the  shock  does  not 
result  in  satisfaction  of  the  Chapman-Jouguet  condition. 

The  structure  of  shock  waves  involving  chemical  reaction  can  in 

0 

some  cases  be  idealized  so  as  to  separate  the  effects  of  viscosity  and 
heat  conductivity  from  those  of  chemical  reaction.  This  model,  which 
is  applicable  if  the  chemical  reactions  are  relatively  slow,  results  in 
a  shock  wave  consisting  of  two  regions;  in  an  initial  region  the  structure 
would  be  described  by  the  equations  employed  above,  perhaps  generalized 
for  variable  specific  heat  In  this  first  region  the  creation  and  des¬ 
truction  of  chemical  species  would  be  zero  so  that  the  clicmical  com- 
I'.c  iGoii  would  be  identical  with  that  far  upstream  of  (he  wave.-...  In  the 
•Second  region,  reaction  would  take  place;  to  some  approximation  it  is 
possible  to  neglect  therein  molecular  transport,  i.e.,  to  neglect  vis¬ 
cosity,  conductivity  and  diffusion  since  lor  extended  reaction  zones  the 
gradients  may  be  sufficiently  small  lobe  inessential.  The  .state  of  the 
gas  at  the  down.stream  edge  of  tlie  first  region  is  chemically  a  non- 
equilibrium  one  witb  the  approach  to  equilibrium  occurring  in  the 
second  region.  This  behavior  has  recently  found  importance  in  ultra- 
velocity  air  flows,  i  e,  .ihosi-  witli  velocities  of  35,000  fl/sec  and  higher. 
In  Ib-is  velocity  range  radiation  from  the  gas  provickus  nn  i rnporiani, a nd 
a.s  the  velocity  increases  the  domi nanliinode  for  energy  ticansfer  to  the 
body:  .since  tiiis  r.'idiaiive  transfer  is  highly  temperature  dependent. 

('•!  .‘'1'*^)'  i.'^.  found  tlias  non-equilibrium  behavu  or  which  occurs  behind 

tiu’  sho;,!;  ln»j  which  n;  .use  icadK  to  rtujic  tvui.p.-ralu  -  er  highc  :•  .than 


equilibrium,  plays  a  dominant  role  in  determining  this  transfer.  Rose 
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and  Teare  provide  a  recent  review  of  tliis  phenomenon. 

For  the  purpose  of  describing  in  general  the  structure  of  shock 
waves  with  chemical  reaction,  the  basic  equations  will  be  presented 
belo\v.  These  equations  should  be  considered  extensions  of  the  Navier- 
Stokes  equations  used  above.  Their  validity  for  either  detonation  waves 
or  for  strong  shock  waves  is  subject  to  at  least  the  same  doubt  as 
discussed  above  but  their  application  is  justified  by  their  simplicity, 
flexibility  and  in  general  by  the  ignoraaice  associated  with  the  details  of 
chemical  kinetic  and  of  molecular  transport  effects.  The  equations 
presented  below  are  at  least  consistent  if  inaccurate;  there  is  Jio  experi¬ 
mental  evidence  known  to  the  authors  concerning  their  validity.  In 
connection  witJi  the  vaJidity  of  the  Navic.r-Stokes  equations  and  their 
extensions  to  detonation  waves  and  strong  shock  waves,  attention  is 
called  to  the  importance  of  the  temperature  dependence  of  the  viscosity 
coefficient  on  the  shock  wave  structure  as  predicted  by  these  equations 
for  the  classical  case  of  homogeneous  gases  with  constant  coefficients 
of  specific  heat  and  =;  3/4,  In  shock  waves  involving  chemical 
reaction  the  temperatures  arc  sufficiently  high  so  that  large  changes 
hi  the  transport  properties  of  the  gas  occur  and  thus  so  that  relative 

The  discussion  of  the  basic  cqtLations  presented  here  is  based  on  the 
standard  treatineuts  of  chemically  reacting  gas  flows,  e,g.,  in 
references  53  and  54 
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thickening  of  the  wave  occurs. 

The  equations  of  conservation  of  mass  and  momentum  (Eqs.  (1) 
and  (2))  are  unaltered  if  reactions  take  place  and  if  it  is  understood  that 
the  viscosity  coefficient  is  now  dependent  upon  composition  as  well  as 
temperature.  However,  the  energy  equation,  which  can  be  written  in 
several  ways,  is  considered  In  the  form; 

N 

ou[(uV2) +h]' -  (kT')' -  (4/3)(muu0' +  (  D.V.h.)  '  =  0  (49) 

i  =  l 


where  the  enthalpy  h  is  now  the  enthalpy  per  unit  mass  of  mixture; 

N 

h  =  )  Y.h,;  p.V,  is  the  mass  flux  of  species  i  due  to  dilfusion  per  unit 
i  =  l 

area  per  unit  time;  h^  is  the  enthalpy  of  species  i  per  unit  mass  including 
the  chemical  enthalpy;  and  where  N  species  are  assumed  present.  Eq.  (49) 
can  be  integrated  once  to  yield  after  .some  rearrangement 


[(u®/2)  Ml]  -[k/(m  c  )]  |[(uV2)i-h]'  -  [  1  -  (4Pr/3)]uu  ' 


t 


y  |(h./m  P.r)M  Y.' [l-KPr  o-V./u  Y.')]|  -  C, 

4-1  1  i  111)' 

1  =  1 


(50) 


N 

\  ■  . 

\vlierc  c  =  )  Y.c  ..  c  .  =  dh./d'r,  and  where  the  Prandil  number  involves 
p  L  i  p,i'  P, 1  1 

i  =  i 

the  niixiure  quantities,  uc^^/k.  Eq,  (50)  indicates  lliai  si,  a-s  previously, 
Pr-3/4.  and  if  PrOjV./i.iY/  =- i  for  all  i,  tlv  prcviou.s,  energy  equation 
with  the  enthalpy  properly  interpreted  applies  to  this  more  g’dicral 
case  of  reactinu  e'Iqw.s  Mow  the 


<■>3 


As  an  c-Xamplc  of  the  temperature  dependence  of  the  !.ranBport  prop- 
;ics  of  air  it  may  be  noted  that  bvCantie  of  ihuj  OnAel  of  iJartial  iOnir.ai 

for  T  <  800CM^K  assd  ak 


the  condiu  livily  of  asr  bcliaves  rbu^hty  as  k"T 
k:  ~  T'’-;’ fo  r  T>BOOO®K.  This  rc.«ults  in  a  sha  rp  dec  rt-n  sc  in  P 


and  s'uuucsi.s  t'nat  some  of  tiie  efiM  ts  refateci  to 


p. 


for  T>8000-'K 
0  described  above  may 


a  r  1  .s  e 


diffusional  mass  fluxes  in  general  arise  from  gradients  of  pressure,  tern 
perature  and  composition  (of.  references  53and.54);n  aeronautical  appUca 
tions  the  diffusional  fluxes  due  to  composition  gradients  are  generally 
assumed  to  be  dominant:  accordingly,  if  a  single  diffusion  coefficient  D 
can  be  assumed,  then 


p.V.  =  -  r.D  Y/ 
11  1 


(51) 


If  pDPr/pi  -Le-l,  then  the  actual  reduction  of  Eq,  (50)  is  realized. 

The  parameter  pD/u  is  the  reciprocal  of  the  Schmidt  number,  Sc,  while 
Pr/Sc  IS  the  Lewis-Semenov  number.  The  approximations  attendant 
upon  Le-  1  are  widely  used  because  of  the  resulting  simplifications  in 
the  gas  dynamics  of  reacting  flows  although  it  is  recognized  that  the 
actual  Lewis  number  in  mixtures  involving  species  of  widely  different 
molecular  weights  (e.  g.  ,  mixture.s  of  atomic  and  molecular  hydrogen 
in  ail')  may  be  greatly  different  from  unity. 


In  the  more  general  case  wherein  realistic  tran.sport  des- 
ciiptiftns  arc  considered,  the  Prandtl  number  can  be  taken  to  be  depen¬ 
dent  on  composition  and  temperature  and  the  diffusional  mass  flu.xes  to 
be  depcndcMU  on  gradients  in  the  flow  variables  in  a  more  accurate 
manner.  In  particular,  the  contribution  to  the  diffusional  mass  fluxo-S 
bs  composition  gradients  are  described  in  matrix  form  as 


In  boundary 
si  mpl duration 


iaver  flows  (ho  values  of  Pr=Sc=:Lc=!  lead 
in  the  energy  equation. 


to  essential 


a. , 
ij 


1  1 


N 


-pD. 


b.  . 
ij 


0  0 . 0 


(5  2) 


Y' 

N 


where  D  is  one  binary  diffusion  coefficient,  i.  e.  ,  for  the  diffusion  of 
r 

one  species  i  in  the  other  j  and  is  used  for  non-dinnensionali'/.ation;  where 

the  a,  coefficients  are  functions  of  the  composition. of  the  molecular 

weights, and  of  the  ratio  of  binary  diffusion  coefficients  to  D^,;wn  mxv  the 

b..  coefficients  are  functions  of  composition  and  molecular  weights;  and 

where  the  last  equation  assures  that  ^ 

i  =  l 

be  seen  that  according  to  this  description  of  the  diffusion  due  to  species 

concentration  gradients^Lhe  flux  of  one  species,  c.  g.  ,  PjV^, depends  on 

all  gradients  Y;,  .  .  .  Y^'  and  not  ^^.Y  /  alone.  Of  cour.so  after  these- 

equations  are  solved  in  the  form  p.V.  ^  c^.Y/,  a  mixture  diffusion 

ml 

cocrficioiil  n  ,  can  be  introduced  so  that  t'',;  V;--- P  but  there  i.s 

only  a  tonuat  advant:u'c  in  dbing  so. 

It  IS  perhaps  worth  noting  at  this  point  that  in  conjunction  with 
the  energy  equation  some  statement  concerning  llic  thermodynamic  slate 
of  the  species  in  the  flow-  must  be  made  The  analysis  of  Detho  and 
Teller’’'  showed  in  a  homogeneous  mixture  ofdiato  -’T  gasc.s  that  the 
iraivsiational  degrees  of  freedom  were  quickly  coihlibraud,  the  rotational 
were  equilibr-ated  next  and  also  rriaiively  quickly,,  but  that  the  vibratiouul 


degrees  of  freedom  were  excited  only  after  many  collisions.  As  a 
reasonable  approximation  for  gas  mixtures  it  is  customary  to  assume 
thermodynamic  equilibrium  with  respect  to  all  degrees  of  freedom 
except  the  vibrational  ones,  ^  For  these,  rate  laws  yielding  .  in  terms 
of  slate  variables  are  required.  The  inclusion  of  this  latter  sophistication 
is  vitiated  in  chemically  reacting  flows;  the  effect  of  coupling  between  the 
chemical  effectiveness  of  collisions  and  the  degree  of  vibrational  or 
rotational  excitation  is  so  complex,  and  available  estimates  of  reaction 
rates  in  practically  interesting  gas  mixtures  so  poor,  (hat  it  is  generally 
not  consistent  from  an  accuracy  point  of  view  to  describe  non-equilibrium 
behavior  in  internal  degrees  of  freedom.  In  the  further  discussion  it  will 
be  assumed  that  full  thermodynamic  eqtiilibrium  prcwails.  Thus  tlie 
forward  and  reverse  reaction  rales  will  be  related  by  the  equilibrium 
constant  in  each  reaction  step. 

The  spccie.s  conservation  equations  arc 

duY/  =  -(o.V.)'  +  -iv.  i=l,?..,.N-l  (53) 

where  w.  is  the  mass  volumetric  rate  of  production  of  apecie.s  i,  i.e.  , 
gms/cc  sec.  Note  that  only  N-1  such  equations  must  be  conlSidei'cd  .since 

Hansen  el  al  have  recently  .suggesled  that  at  high  tem pc,*  ralu  re s 
rotational  relaxation  m.iy  be  si  enifica nt,  Talbot  and  Scala^  ^  at  the  same 
time  studied  iheo.reiically  .shock  structure  in  a  di-atomic  gas  with 
rotational  relaxation:  tlteir  main  motivation  appears  lo  be  the  d,ete>r- 
minatiqn  of  a  suitable  inoclel  ,  i.e,  bulk  visco.sity  or  relaxation,  for 
lagging  internal  degree's  of  freedom. 


Eq,  (1)  assumes  overall  energy  conservation.  Eq.  (53)  indicates  that 
the  gradient  of  concentration  of  a  given  species  at  a  generic  point  in  the 
flow  depends  on  diffusion  and  on  chemical  reaction.  The  creation  term 
w.  is  determined  by  the  mechanism  of  the  chemical  reactions  involved; 


for  a  system  involving  L  reaction  steps 


_  l  1  . 

w,  =  W.  )  (vf.-v/.)k,  II  (Y./W  )\i  G  p  ^ 
1  i  G  '  ij  ij  j  1  1  g  J 


where  v"  and  v'  are  the  stoichiometric  coefficients;  k  is  the  specific 
ij  iJ 

rate  constant;  G.=  0  is  the  equilibrium  condition,  i.e,, 

J 


G.  =  1  -  (p  Vk  J  I!  (Y./W  ) 

j  c ,  j  .  _ ,  I  i 


//  / 
V.  .  -  V.  . 

\  IJ  bb 


n  -  \  v'  ;  m  -  )  v"  -  V.'.  ;  K  .is  the  equilibrium  constant  for  reaction 
j  “  G  ij'  j  -G  ij  ij  c,  J 

i=l  i=l 

step  i  in  terms  of  molal  concentrations,  and  where  the  j  inaction  step 


is  syrnbollically  given  by 


k.  N 


)  v.'.M. — V.'".  M. 
Zj  1  )  1  i_  1  )  1 


Eq.  (53)  can  be  integrated  only  formally  because  . of  the  creation 


te  rm,  i .  c.  , 


Y.  I  (o.V./m)  =  ;  {\v. /m)  clx  t  constant  (56) 

i  i  i  ^  ‘ 


'rhi.s  eouation  indicates  al.so  the  influence  of  diffosion  ma.ss  flux' and  of 


chemical  reaction  on  the  distribution  of  mass  fractions. 

The  analysis  of  composition  is  simplified,  if  the  eqviations  of 
element  conservation  are  considered.  Suppose  in  the  general  system 
involving  N  species  there  are  M  elements;  the  conservation  of  elements 
requires  that  throughout  the  flow 

N 

L  =  0  j  =  1.2,  ...M  (57) 

i=l 


where  p,.  is  the  number  of  atoms  (or  molecules)  of  clement  j  in  one 
molecule  of  species  i.  .Eq.  (5  7)  suggests  an  element  mass  fraction 
Y.,  that  is,  the  mass  density  of  element  j  in  all  species  per  ma.ss  of 
mixtuJ-c,  given  by 
N 

Yj  ,.,M  (58) 

i  =  l 

1'  rom  .Eq.  (53)  with  Eqs.  (57)  and  (58)  considered,  tlie.re  is  obtained 


N 


i  =  ! 


(59) 


or  upon  integration 


N 

raY.  f  )  u.  ,W.r;.V,.AY.  =C 

J  L  ij  j  1  /  1 

1  =  1- 


I  he  set  of  equanons  providing  the  formal  description  of  shock 


(61) 


structure  in  reacting  gases  is  completed  by  an  equation  of  state 

N 

P  =  pR  T  y  Y./W. 

0  r  1. 

i  =  l 

where  is  the  universal  gas  constant;  by  equations  yielding  the  co¬ 
efficients  of  viscosity,  specific  heat  and  heat  conductivity  as  functions 
of  composition  and  state;  and  by  descriptions  of  the  functional  dependence 
of  the  reaction  rates  and  equilibrium  constants  on  temperature. 

The  entropy  distribution  through  the  wave  may  be  computed  oiice 
the  distributions  of  composition  r.nl  .o')  i  r  o  vl  c  va  riablcs  are  kjiown;  the 
entropy  per  unit  mass  of  mixture  is 

N 

S  r  -iWjW')  (n{p\V)  I  ^  hUY./W^I]  |u2| 

i=l 

is  tho  stamiard  ouiropv  of  hin-tics  i  per  unij^uiass  of  .sp.  *  ic  s 

i  and  •  li  r.- W  i-  Jill  inixi'ire  n-.olfC'ilar  weii^l-.t.  W‘  -  (Y  AV  > 

I  I 

1^1 

wunsioei  now  the  application  ol  these  conservation  cquatiorns  a( 
poinl.s  ol  eMilonmiv  in  the  flcnv,  i  e.  ,  at  x  =  ±Coi  lltr  in;,';.-  and  (notroTtun. 
conservation  as  before  lead  to 

0  U  ;;;  m  -  u 
it  3  ? 

mu  -f  p  =  rnu  f-  p 

*  C  i'. 

and  (';0)  VI'  Id  r.  sy.<  .  tiv'dy 

=^2 


wJlilr  ncs,  (aO), 


(63a) 


3 

(u  jZ)  +  h 
1  ^  1 


(u  /2)  -h  h 
z  z 


/  /  n  ■<  » 

(O  JD) 

(63c) 


since  vanishes  for  all  i  at  x  ~  :!  co  (cf,  ,e,g.  ,  equation  (52)). 

If  the  initial  aird  final  states  correspond  to  cheinica.1  as  well  as  thermo¬ 
dynamic  equilibrium,  then  the  equations  of  mass,  momentimi,  energy 
and  element  consei'vation  (Eq.  (63c)),  i.e,,  those  yielding  simple 
algebraic  relations  relating  the  flow  and  state  variables  at  1  luid  2, 
can  be  supplemented  by  ^u^  equation  of  state  and  by  a  sulficient  number 
of  the  equilibrium  conditions  G.^^O;  c.g.  ,  for  a  chemical  system  with 
N  Species  and  M  elements  N-M  equilibrium  conditions  must  be  specified. 
Thus,  if  all  conditions  corresponding  to  1  arc  known,  the.  aforementioned 
equations  determine  all  conditions  at  2;  in  the  case  of  detonation  waves 
satisfying  the  C-.T  condition,  the  velocity  u  is  also  unknowi  but  the  re¬ 
quirement  that  u  equal  the  Speed  of  soinrd  at  plus  infinity  completes  lire 
system  of  equations.  Eihropy  considcratiojis  may  have  to  be  invoked 
in  order  to  select  the  practical  dowTistream  co7Kiitions  if  multiple  j)os- 
sibilitics  therefor  arise. 

In  considering  the  shock  wave  Structure  it  can  thus  be  assumed 
that  the  end  poir.ts,  i.e,  ,  conditions  at  1  and  2,  arc  knov/n.  The  range 

2  a 

Sec  refc rence  4B  for  a  discussion  of  the  app ropriatK  speed  of  sound  in 
view  of  the  existence.,  of  the.  two  limiiitjg  ca.Hcs  of  "fronen’-'  and 
*' equilibrium  ’  speeds  of  sound. 
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of  independent  variable  can  be  changed  from  ±co  to,  e,  g,  ,  0  to  1  and  the 
numerical  analysis  can  be  simpliiied  by  introducing 


<  =  (u  -u]/(u  -u  ) 


Then  the  diffeiential  dx  is  eliminated  from  the  energy,  species  and 
element  conservation  equations  since  '=d/dx=  -  u '(u  -u^ )' ^  d/di  and 
since  ti  '  can  be  eliminated  by  the  momentum  eq-;ation  integrated  once. 

This  transformation  of  the  independent  variable  is  singulaj-  at  =  0,  1  so 
that  the  asymptotic  behavior  of  the  solutions  for  the  dependent  variables 
must  be  developed  at  these  points.  After  the  solutions  in  terms  of  I 
are  obtained,  the  monentum  equation  can  be  integrated  a  second  time 
by  quadrature  to  obtain  x  “x(T)  subject  to  same  condition  fixing  the 
origin  of  the  x  coordinate. 

To  give  some  indication  oJ  the  effect  of  cli emica I  reaction  on  the 
structure  oi  a  sliock  vravc  consider  Fig.  10  which  i.s  taken  with  modifienljon 
and  extension  from  reference  d?."'  Tiie  conditions  far  upstream  and  down¬ 
stream  of  the  wave  are  shown,  Characteristic  of  sliock  structure  involving 
finite  rate  chemistry  is  the  aforemctuioned  overshoot  of  tlie  temperature  so 
that  within  the  wave  temperatures  exceed  1/,  The  flow'vanablos  re!aicci  by 
the  momentum  equation  on  the  other  hand,  i.  e,  ,  velocity  and  pressure 
behave  much  as  in  the  classical  shock  wave  Shown  in  Fig.  IQ  is  ihe 
distribution  of  the  total  mass  fraction  of  di  s  soc  lated  species,  m  this 

In  it*i(‘i  cnce  •!  i  (jaila tae.s  and  Bloom  .eon.side  red  a  mixiure  orudeal 
dissociating  gases  behaving  ii]!;c  oxygen  and  nitrogen,  a .s smned  Pr=5/-i 
and  tor  all  species,  neglected  diffusion  in  the  species  consm-vaiiun 

equation;  and  assumed  thermodynamir  equilibrium,  Jniegratmn  was  ^ 
started  at  conditions  corresponding  to  -ai,  but  at  a  finii,-.  valup  of  x 


o 


case  0  and  N;  clearly,  an  analysis  based  on  the  von  Neumann-Doring- 
Zeldovitch  model  (cf.  reference  48)  of  a  detonation  point,  i.  e.  ,  on  a 
shock  wave  without  chemical  reaction  followed  in  the  downstream 
direction  by  a  region  of  chemical  reaction  appears  to  be  valid  according 
to  the  assumptions  of  this  analysis.  More  realistic  descriptions  of  the 
transport  phenomena  may  alter  this  conclusion. 

The  entropy  distribution  in  the  wave  has  been  computed  from 
Ecj.  (62)  enrploying  the  temperature,  pressure  and  composition  piofiles 
given  by  reference  47  and  is  shown  in  Fig.  10  in  terms  of  (S-S^)/R^. 

It  will  be  noted  that  the  relatively  simple  entropy  distributions  arising 
in  classical  shock  wave  structure  no  longer  arise;  this  is  undoubtedly  due 
to  the  interplay  between  heat  conduction  and  chemical  kinetics  in  this  more 
general  case.  In  the  initial  portion  of  the  wave  structtire  leading  to  a 
plateau  in  the  entropy  the  effect  of  heat  conduction  predominates  so  tliat 
the  simple  behavior  of  the  classical  shock  wave  begins  to  appear.  How¬ 
ever,  in  the  more  downstream  portions  of  the  wave,  chemical  reaction, 
i.c.  ,  dissociation,  takes  place  and  leads  to  further  increa.ses  in  entropy. 

A  slight  maximum  in  the  entropy  occurs  at  ?  -  4,  roughly,  but  no  con¬ 
clusion  as  to  the  significance  or  generality  of  this  result  can  at  j^resent 
be  d  rawn .  ^ 


After  coinpletion  of  ihi.s  manuscript,  a  sin-vey  of  the  problem.^ 
connected  by  shock  strr.cture  by  Talbot^  ^  appeared. 
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